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ABSTRACT 



Uniform and isotropic mathematical models of the expand- 
ing universe usually predict an initial singularity of infi- 
nite mass density and space curvature. To study possible 
mechanisms which v/ould avoid the occurance of these singular- 
ities, non-uniform cosmological models based on Einstein's 
field equations are investigated in v;hlch random perturba- 
tions of long wave lengths are superimposed on the Robinson- 
Walker metric of the unperturbed models. Techniques of fluid 
turbulence theory, used to describe random fields by a hier- 
archy of central moments of the random perturbations, are 
applied to describe the dynamics of these moments. For the 
case of small perturbations the hierarchy is truncated and 
solutions are found. The solutions are either grov'ing or de- 
caying perturbations leading to r'^ extra terms in the usual 
cosmological equations for the curvature radius R. The re- 
sult agrees with the small perturbation Fourier series expan- 
sion analysis which exists in the literature. Based on the 
upper limit of the anisotropy of the 3° K background radia- 
tion, the growing perturbation model predicts a maximum 
expansion even for k=0. Euclidean spaces. The decaying per- 
turbation solutions give extra terms of the form 1/R^ with 
m>4 in the cosmological equations and indicate that the 
mechanism of long v/ave random perturbations may prevent the 
original singularity and make oscillatory models possible. 
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I . INTRODUCTION 



Perhaps one of the most exciting discoveries of the 

Twentieth Century thus far is Hubble's now famous discovery 

( 1936 ) that the relative shift In wavelength — in the 

A 

spectral emissions of each star is proportional to its 
astronomical distance L [1]. Modern cosmology generally 
recognizes this phenomenon to be the consequence of an 
"expanding universe" whose dynamics are determined by the 
Einstein field equations [2, 21], The concept of an expand- 
ing universe presents various problems to those v.'ho are 
involved with the systematic study of our universe as an 
entity. Because the term "expanding universe" implies that 
all distances between all points are on the average increas- 
ing, it is natural to ask if our universe was ever in a 
condition in which all lengths were zero. It is also natural 
to ask if this so-called expansion is homogeneous and iso- 
tropic at all points in the universe. 

If the ansv/er to the first question is yes, then it is 
apparent that sometime in the distant past the mass density 
of the universe was infinite. But V/heeler [20] has pointed 
out that such a state is rather dubious. He calls attention 
to the fact that it is possible to construct a "length" L* 
by a combination of the gravitational constant G, the 
light velocity c, and the Planck constant ti such that 
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1.6 X 10 



cm. 



( 1 - 1 ) 



L* 






-33 



He suggests that this length might be a characteristic 
length that is associated with a quantum theory of gravita- 
tion. He points out that because of the smallness of this 
length, quantum gravitational effects, though negligible 
now, could in principle be important during the very earli- 
est stages of the expansion of our universe [5]. Hence, 
it is believed that the undesirable singularity associated 
with the infinite density might be avoided. One reason, 
he says, why such quantum effects might be significant 
during the highly contracted phase of our universe is that 
the estimated cross section for pairs of gravitons to produce 
pairs of electrons and pairs of photons is comparable when 
the wavelength of the compacted gravitational radiation, 
which is assumed to be present, is of the order of L* . In 
order to describe such a situation, it is clear that a 
quantum theory must be used to treat this obviously non- 
classical phenomenon [20]. 

Various workers have begun the task of incorporating 
quantum effects into gravitational theory. Some investiga- 
tors (Misner [l6] and Parker [1?] for example) do this by 
various canonical procedures, but others have taken a differ- 
ent approach. V/oehler [3^] suggests that it might be fruit- 
ful to postpone a correct quantum theory of gravitation and 
try a semlclasslcal approach. A way he suggests to do this 
is to presume that the field variables g^^ (called the 
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metric tensor) of Einstein's field equations randomly fluc- 
tuate with position and are characterized by a "mixing 
length" L, v;hlch for quantum fluctuations v;ould take the 
value L* • By considering the average of the field equa- 
tions, he suggests that it might be possible to avoid singu- 
lar solutions if the extra source term, v;hich appears as a 
consequence of the averaging and the nonlinearity of the 
field equations, is dominant during the early phase of high 
contraction of our universe. Ginzburg and other co-v/orkers 
[ 7 ] use a similar approach. They attempt to estimate the 
functional form of this extra source term by expanding 
this source term in a series expansion of the average curva- 
ture tensor. It is interesting to note that both the tech- 
niques of Woehler and Ginzburg are reminiscent of an earlier 
technique used by Welton [33] to estimate the Lamb Shift 
by a semlclasslcal model that also involves fluctuations, 
but of the position variable of the electron. 

However, the point in the above discussion that is 
germane to the introduction of this thesis is that there 
are today serious considerations of using a randomly fluctu- 
ating metric tensor g^^ as a prelude to developing a realis- 
tic quantum theory of gravitation. 

There is also another aspect of cosmology that requires 
the use of a nonuniform fluctuating metric tensor 
This consideration is brought out by the second question 
that asks if the expansion is everyvjhere homogeneous and 
isotropic. It is a matter of history that one of the first 
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attempts to describe the universe from a cosmological point 
of view assumes that all matter and all space are homogene- 
ous and isotropic on a sufficiently large scale. This 
theoretical assertion is called the Cosmological Principle 
and v;as formulated in 1935 [2, 21, 22]. 

Although the use of the . Cosmological Principle appears 
regularly in the literature, its criticism is still very 
much alive today. De Vaucouleurs [31], in a recent paper, 
points out that because current observations have as yet 
failed to demonstrate the existence of that sufficiently 
large scale beyond which space is homogeneous and isotropic, 
it seems unlikely that it v/ould be easily found (if such a 
length scale exists at all) by a slightly deeper astronomical 
survey . 

Because of this and other arguments, many workers have 

been led to consider nonuniform cosmological models. One 

» 

way this is done is to consider the dynamics of an anisotrop- 
ic universe which has tv;o equivalent and one inequivalent 
directions in space. The dynamics of such models can be 
determined exactly by solving the appropriate field equations 
[ 27 ]. Another method to construct nonuniform cosmological 
models is to consider small linear perturbations away from a 
spatially homogeneous and Isotropic cosmological metric . 

Sachs and VJolfe [24] do this and are able to express these 
perturbations as a function of time. Their solutions contain 
some arbitrary functions of the three spatial coordinates 
which are determined by initial conditions of a given 
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arbitrary perturbation field. These functions then are 
generating functions which correspond to a given nonunlform 
model . 

Again the Important point for the reader of this paper 
to note is that cosmological models v;lth a variety of non- 
uniformities are being actively considered by modern cosmol- 
ogists. These nonuniformities may be zero point quantum 
mechanical fluctuations or they may be fluctuations due 
to a large scale uneven mass distribution. In fact, Ginz- 
burg and others [7] point out that both types of fluctua- 
tions are described by the same formalism, a metric tensor 
that fluctuates randomly about its average value. Matzner 
[l4] points out that this formalism is appropriate for a 
universe whose average background is uniform but contains 
gravitational radiation of an arbitrary size L. Clearly, 
some techniques are needed v;hlch describe these random per- 
turbations of arbitrary magnitude and size L in some natural 
way . 
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II. MOTIVATING PHILOSOPHY OF THE MODEL 



In 193^5 Milne and McCrea chose to describe our universe 
In terms of a Newtonian model. This model treats the uni- 
verse as a streaming, uniform fluid. Surprisingly, they 
showed that such a model Is In many respects equivalent to 
the models of relativistic cosmology [2]. Their model has 
the distinct advantage of being mathematically much less 
complex than models developed within the framev/ork of gener- 
al relativity. Another advantage Is that fluid flov/s have 
been well studied and It v/ould seem, therefore, that some of 
the fruits of fluid mechanics might be transferable to the 
field of cosmology. 

Such a possibility has. In part, been the motivation 
for preparing this dissertation. In particular, the study 
of turbulent fluid flov;s which deals with random perturba- 
tions of fluid flow fields would seem, on the basis of the 
above discussion, to be relevant to the study of random 
perturbations of the "flow field" associated v;lth the 
universe . 

It Is well known that fluid flov;s obey the Navler- 
Stokes equation. This equation relates the velocity field 

-h 

V to a pressure field P associated with the flow. It Is 
also v/ell known that If both the velocity and pressure 
fields are decomposed Into an average field plus a fluctua- 
tion field such that 
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(2-1) 



V = <V> + V 
P = <P> + p 



where < > denotes "average," then the Navier-Stokes equation 
can be placed into a form v;hich is symbolically identical 



with the Navier-Stokes equation for V and P but with V 
replaced by <V> and P replaced by <P>. Owing to the non- 
linearity of the Navier-Stokes equation, additional source 
terms v/hich involve the correlations <v^v^> appear (v;here 
i and j denote the i and components of the fluctuation 

v) [^]. These extra terms are called the Reynold Stresses 

and may be Interpreted as additional stresses which act 
-> 

upon <V> . 

There are a variety of methods that fluid dynamlcists 
use to determine the functional form of the Reynold Stresses 
An important way is to generate a hierarchy of moment equa- 



tions in V by multiplying the unaveraged Navier-Stokes 
equation with all orders of V and averaging. Because of 
the nonlinearity, each new equation so generated contains 
a moment of one order higher than v;hat would be present 
if the Navier-Stokes equation were linear. Although any 
number of moment equations could in principle be generated, 
there always exists one more unknown than there are equations. 
This problem is called the closure problem of turbulence 
theory. Much research has gone into finding ways to replace 
this indeterminate, infinite set of equations by a plausible 
determinate, finite set so that useful information about 



turbulant flow patterns might be obtained [3]- However, in 

li» 



spite of the difficulties associated with the closure prob- 
lem, it should be pointed out that the turbulence approach, 
which is a sufficiently general formalism, can in principle 
describe any sort of random perturbation field of any magni- 
tude. (A short summary of turbulence theory of fluids is 
given in Appendix A.) 

The Einstein Field equations are likewise nonlinear 
differential equations. Their field variables are the 
metric tensor g^^ and the energy-momentum tensor 
Because certain aspects of the field equations of general 
relativity resemble those of fluid mechanics, it is not 
unreasonable to presume that the Einstein field equations 
could be expressed in terms of the average metric tensor 
the average energy-momentum tensor some 

residual terms that appear in a way analogous to the Reynold 
Stresses ov;ing to nonlinearities. It might also be possible 
to generate a moment hierarchy which can be used for the 
determination of these extra terms if an appropriate closure 
technique could be found. 

It is this speculation therefore, based on the similar- 
ities betvieen fluid mechanics and general relativity, that 
serves as the motivation for constructing such a technique 
which might be useful in solving perturbation problems of 
general relativity. The remainder of this dissertation is 
a reply to this conjecture. 



15 



III. INITIAL DESCRIPTION OP THE MODEL 



Given that it is possible to generate a hierarchy of 
meaningful moment equations from the Einstein field equa- 
tions, it is necessary to truncate them in some manner. 

Then a solution for the "Reynold Stresses" can be determined 
which will be found to contain second order moments in the 
perturbed metric and its derivatives. Just as in fluid 
turbulence theory, the closure of the moment hierarchy is 
nontrivial because it requires special restrictive assump- 
tions concerning the nature of the random perturbation 
field and does not just "fall out" of the theoretical fra 
framework . 

In order to see how this truncation might be accomplished, 
consider first the following three important truncation 
procedures v;hlch are used in fluid turbulence theory to 
close its unlimited set of equations v/hich contain moments 
of all orders of the velocity v [3l* 

1. Phenomenological approximations - These involve postu- 
lating (on dimensional grounds or by intuition) an 
extra relation between various moments to bring the 
number of equations up to the number of unknowns. 

2. Higher-order-moment-discard approximation - This tech- 
nique simply discards the higher order moments because 
of smallness or some other reason. This reduces the 
number of unknowns to the number of available equations. 
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3. Cumulant-discard approximation - This procedure relates 
a higher order moment to products of lower order moments 
and thereby effectively establishes an extra equation. 

All these closure methods are used in fluid turbulence 
theory and so all might conceivably be acceptable candidates 
for closing a moment hierarchy in general relativity as well. 

If it is desired to describe cosmological models with 
a ’’small" degree of anisotropy and inhomogeneity, then it 
would seem that the appropriate truncation would be number 
two, the higher-order-moment-discard approximation. This 
choice is reasonable because it v/ould be expected that a 
moment of order n+1 v;ould be negligible compared to a moment 
of order n. Such weakly perturbed models are considered, 
for example, by Sachs and V/olfe [2^] as previously stated. 

However, not all cosmological models are adequately 
described by small perturbations. In a recent paper by 
Woehler [3^], a model is presented which assumes quantum 
fluctuations of the metric tensor to be as large as the 
average metric tensor itself. Certainly the higher-order- 
moment-discard approximation is not adequate to apply to this 
model because all moments of all orders are important, espe- 
cially near the singularity. For this type of model, the 
first or third truncation alternative might be useful. 

Therefore, an attempt to construct a general perturba- 
tion technique based on a moment hierarchy scheme which is 
able to treat fluctuations of all shapes and sizes is dis- 
couraged by the fact that the method of truncation depends 
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upon the type of fluctuation that is involved. In fact, 
even the solutions of the moment hierarchy could, in general, 
be quite different owing to the fact that different moment 
equations could result from different methods of truncation. 

So it seems that a reasonable viay to proceed is to 
consider a model with a certain "kind" of perturbation 
field and then form and truncate its moment hierarchy by a 
technique that is compatible with the type of fluctuations 
assumed. Thus, a two fold purpose would be achieved. First, 
a model would be constructed which describes a special non- 
uniforrnity about which conclusions can be drawn that may be 
of value to cosmology. Second, and perhaps more signifi- 
cant, would be the demonstration that this procedure does 
work and is therefore relevant to perturbation problems 
of general relativity. 

In order to best serve the process of demonstrating 
the relevancy of the moment hierarchy technique to general 
relativity, a simple nonuniform model is used. It is simple 
in the sense that only one function is required to describe 
the perturbation field rather than ten which are needed to 
treat a general perturbation problem. This random function 
is also assumed to contain random perturbations which are 
small compared to the average field. Thus, the previously 
discussed higher-order-moment-discard approximation will 
be used for the truncation of the moment hierarchy which 
will be constructed in a v:ay analogous to fluid turbulence 
theory. Even for this relatively simple model, the 
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mathematical manipulations are quite extensive. It v/ill be 
shown that this particular model is relevant to the study 
of long wavelength perturbations and deformations of the 
size which approaches the length characteristic of the 
universe itself. Thus, the model, though highly restrictive, 
has a special relevancy and direct applicability to 
cosmology . 
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IV. METRIC TENSOR 



A. SPECIAL FORM FOR THE MODEL 

The process of building the model begins v;ith the formu- 
lation of the metric tensor This can be accomplished 

by giving first the equation for the square of the line 
element ds relative to a given coordinate system and then 
computing by the relation 

ds^ = g dx'^dx^ . (^-1) 

To this end, consider first the line element for uniform 
cosmological models • It is called the Robert son-Walker line 
element and in spherical coordinates is given by [1] 



where 



and 



ds^ = (dx®) - e*^^^ 



da^ = dr^ + r^de^ + r^sin^6d({)^ 



G(x®,r) _ ^ 



R " (x®) 



1 + 



k r‘ 



^ “7 

r „ 



(^- 2 ) 



( 4 - 3 ) 



( 4 - 4 ) 



The parameter is an arbitrary normalization constant 
with dimensions of length for the radial coordinate r. 
k is the curvature characteristic which may take only the 
values -1, 0, and +1. The time coordinate is x® and the 
three spatial coordinates are r, e, and 4). The function 
R^(x°) is called the scale length. 

u 
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Adler [1] points out that the spatial part of the line 
element of Eq. (^-2) describes a three dimensional geometry 
which is homogeneous and isotropic with a uniform radius of 
curvature R^. Associated with the three values of the 
curvature characteristic k are three intrinsically differ- 
ent kinds of uniform geometries. The first, for v;hich k 
is equal to zero, is the familiar Euclidean geometry. This 
geometry is based upon the postulate that through any point 
outside of a straight line, one, and only one parallel line 
can be dravm. The case v/here k is equal to +1 corresponds 
to the uniform geometry first discussed by Riemann vjhich 
is based on the postulate that any two straight lines in a 
plane alv/ays intersect each other. Finally, the case vjhere 
k is equal to -1 corresponds to the noneuclidean geometry 
of Bolyai-Lobachevaki . This uniform geometry is based 
upon the postulate that through a point outside of a straight 
line, an unlimited number of straight lines can be dravjn 
which do not intersect the given line and which lie between 
two straight lines that do intersect the given line at 
infinity [26]. 

Now, if the function GCx^r) were replaced by an arbi- 
trary function M(x ° ,r , 0 , <j) ) , the resulting line element would 
no longer be necessarily uniform. In fact, a positive 
function Rp^(x® ,r , 0 ,<{) ) could be defined such that 



^M(x« ,r, 0 , 4 ,) _ 



k r2 

^ ^ J — 



(^- 5 ) 
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Thus, at every point x^, a tangent Robertson-V/alker geometry 
could be constructed v;ith a radius of curvature equal to 

U 

the function R.,(x^) . Hence, the function R„ could be viev/ed 
in M 

in this sense as a "local radius of curvature" or "local 
scale length." 

As stated previously, one of the objectives of this 
dissertation is to demonstrate the applicability of a prin- 
ciple of fluid turbulence to perturbation problems of general 
relativity by means of a simple model. Clearly, a model 
whose geometry is described by the line element of Eq . (4-2) 
with G replaced by M is simple because the entire nonuni- 
formity is described by the single function M (in general, 
ten independent functions are needed to describe a nonuni- 
form metric Sy,^(x^) because g^^ is symmetric with respect 
to u,v). Thus, such a model is a convenient choice for 
such a demonstration. It is therefore assumed that the 
line element of the model is given by 

ds^ = (dx°)^ - . (4-6) 

Because this model is to be described by statistical 
methods, the value of the function M at any given point 
must possess a certain amount of randomness. That is, it 
must be viewed as a stochastic process or random field. 

Hence, M is a function not only of the four space-time 
coordinates (x°,r,0,<{)) but also of an ensemble parameter C 
which labels each of those functions M that are possible 
for a given experiment. Hence, Eq. (4-6) describes not just 
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one geometry, but an entire family of geometries each of 
which must obey the Einstein field equations. 

The notion of an experiment in the study of stochastic 
processes is well defined. See, for example, Papoulls [l8]. 
However, for the purposes here, there is an interesting 
connection betv;een an experiment (as used in statistical 
theory) and a modern concept called superspace. Wheeler 
[5] defines superspace as the manifold each of v/hose "points" 
represents one three-geometry. A given submanifold of 
superspace is therefore a history of the geometry undergone 
by physical space. Hence, the term experiment as used in 
this paper would indicate an observation of the "state 
vector" of a possible three-geometry in superspace at a 
given time. Certainly, the set of possible state vectors 
for a given experiment is severely limited by Eq. (4-6). 

If the tensor form of the line element given by Eq. 

(4-1) is compared v/ith the line element given by Eq . (4-6) 
with 

xo = x° 

X ^ = r 
x2 = 0 

x3 = (j) , (4-7) 

it is easily found that the nonvanishing elements of the 
metric tensor for the model must be 



S 00 ^ 
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(4-8) 
cent . 



g22 




g33 




i^-Q) 



B. DECOMPOSITION OF THE METRIC TENSOR INTO ITS AVERAGE 
AND FLUCTUATING COMPONENTS. 

In fluid turbulence theory, each of the random field 
variables (pressure and velocity) of the governing Navier- 
Stokes equation is decomposed into its average and fluctuat- 
ing components as described in Chapter II. This same 
decomposition may be performed upon the random metric tensor 
g which is the field variable of general relativity. 

V 

That is, g^^ may be written as 



where <g^^> is the average of g^^ and f^^ is its fluctuating 
component. This average is understood to be an ensemble 



the value of the random field from each member of the 
ensemble of possible states of the system which occurs 
with a certain probability are summed and divided by the 
number of ensemble members. As this number tends towards 
infinity, this sum tends tov/ards the average value. (The 
concept of average is discussed in more detail in Appendix B. ) 

As seen in Eq. (4-8), the randomness of g is contained 
^ ' yv 

M 

in the stochastic process e . This process can also be 
decomposed at each point in space-time into two functions 




(4-9) 



average. That is, to take an average at each point x^ 



such that 



(^ 4 - 10 ) 



M M . ^ 
e = <e > + 

M M 

where <e > is the average of e and is the fluctuation of 



M 

e about its average. By Eq. (^i-8), the nonvanishing com- 
ponents of <g > and f must be 
yv yv 

<Soo > = 1 

M 

<gii > = - <e > 

<g22 > = - r^<e^> 

-M 



<gs 3 > = - r^sin^0<e > 



111 - - 
f22 = - r^f 



M 



f 33 = - r ^sin^ef 



M * 






Another stochastic process t(x^) can be defined as the 

M 



ratio of the fluctuation of e to its average. That is, 
f 



T E 






M 

<e > 



(4-12) 



If f^^ is expressed in terms of x, its nonvanishing elements 



become 



f 11 = - <e >T 
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2 M 
r^<e >x 



M 



fj3 = - r2sin^0<e >x . 



(4-13) 



Because the average of fj^j must be zero according to Eq. (4-10), 



it follows that 



<T > = 0 



(4-14) 
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Also by Eq. (^^-10), it follows that 



T ^ - 1 



(^- 15 ) 



M M 

because e and therefore <e > must be positive. Certainly 

models v^ith small perturbations v;ould be ones for v;hich the 

ratio of the nonvanishing components of f to <g > is 

yv ^yv 

small; that is, for these models, the magnitude of x v/ould 
be very small compared to one. As stated in Chapter III, 
it is this kind of perturbation that v/ill be considered in 
this work. 

1 . Requirement for the Average Component of the Metric 
Tensor 

Suppose it is assumed that the average metric tensor 
<gyv> is Robertson-Walker . Then, the average geometry v;ould 
be homogeneous and isotropic everyv.’here . This assumption, 
which is used extensively in the literature [7, 1^, 2^], 
is convenient because of its relatively simple mathematical 
form. Furthermore, in the limit of zero fluctuations, the 
model becomes determinate and the geometry is automatically 
uniform. 

Thus, the nonvanishing components of the metric tensor 
<gyv> for the model are given by 



<goo > 
<gii > 
<g22 > 
<g33 > 



- e 



G(x° ,r) 






2 . 2 o G(x® ,r) 
- r sin Oe •’ 



(i]-l6) 
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Because the average metric Is assumed to be everyv/here homo- 
geneous and isotropic and is given by the above equations, 

M G 

it follows that <e > = e . By virtue of Eqs. and 

(^-5)j it can be easily seen that 

,r ,0 , 4 ) )> =Rg^(x®) (^-17) 



where Rj^(x® ,r ,6 ,4» ) is the already discussed local scale 

length of the random metric and R^ is the universal scale 

length of the average metric. If AR^ (x ® ,r , 6 ,<J> ) is the 

2 2 

fluctuation of R.. from R_ such that 

rl u 



then by Eq. 

T = 



(^- 12 ), 

AR^ 

2 



it follows that 



(^-19) 



Thus, T can be Interpreted as the ratio of the fluctuation 
of the local radius of curvature squared to the average 
radius of curvature squared. 

2 . Requirement for the Fluctuating Component of the 
Metric Tensor 

In fluid turbulence theory, a convenient require- 
ment concerning the statistics of the perturbation field is 
the requirement of stationarity [ 133 * If this statistical 
quality is borrovjed and applied to the present model, some 
simplifications can be expected by assuming the statistics 
of the turbulent field t to be stationary with respect to 
a change of coordinates on a spatial hypersurface of constant 
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X®. Clearly, the statistics may not be assumed to be 
stationary with respect to the time coordinate x® because 
this v/ould imply a model universe that is not in a state 
of dynamic evolution (i.e. expanding). In Appendix C, the 
implication of this concept is discussed and the important 
formula given below is derived. 

9.<h[T(x^)]> = 0 (iJ-20) 

J 

where j = 1, 2, 3 and h(T) is any function of the stochastic 
process x. Clearly this is a very useful formula for by it 
many bothersome terms are eliminated. 
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V. EINSTEIN FIELD EQUATIONS FOR THE FLUCTUATING METRIC 



A. DESCRIPTION OP THE FIELD EQUATIONS IN GENERAL 

Just as the Navier-Stokes equation relates the velocity 
field to the pressure field in a fluid, so do the Einstein 
field equations relate the geometry of space to the energy 
content of that space [1]. Being a tensor theory, the 
field equations are expressible in a form independent of any 
coordinate system. The geometrical variable of the field 
equations has already been identified as the metric tensor 
g . The energy content variable is called the energy 

p V 

momentum tensor T which is developed and discussed later. 

pv 

It should be stated, however, that the field equations for 
the model Itself are not expressed in covariant form because 
a special coordinate system given by Eq . (^-7) has already 
been imposed. 

Starting v;ith the metric tensor the field equations 

can be constructed as follows. (This procedure is given by 
Adler [1] and is used in the subsequent sections to build 
the field equations specific to the metric tensor of the 
proposed model.) 

First, the metric tensor g is substituted into the 
’ pv 

follov^;ing expression. 

[uv.6] = I - 3gg„^) (5-1) 

where a e — ^ . [pv,6] is the Chrlstoffel symbol of the 

3x^ 
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^11 






first kind and Is a set of 6^1 elements v;hlch are symmetric 
v;ith respect to their y, v indices. 

Next, a quantity , the Christoff el symbol of the 

second kind, which also is a set of 6^ elements, is calcu- 
lated by 

= s“^Cviv,3] (5-2) 



where the identical indices Indicate a summation over all 
cx 

3. r is also symmetric with respect to its lov;er indices. 

yA v 

CX 6 

The object g , called the second rank contravarlant metric 



tensor, is related to the metric tensor g by 

* °p V 

a3 

8 Spe = % • 



(5-3) 



Next, the Ricci tensor R is formed by 

P V 



R = 3 - a 

pv V pa p pv pv pa pa pv 



(5-^) 



It is a set of l6 elements. The Ricci tensor is then used 
to calculate the mixed Ricci tensor R^^ given by 



V 



X6^ 
g H 



6v 



(5-5) 



Next, the Ricci scalar R is formed by the contraction 
of the mixed Ricci tensor R^^ given by 



R = R 



(5-6) 



Finally, the mixed Ricci tensor of Eq. (5-5) and the 
Ricci scalar of Eq. (5-6) are used to form the mixed Einstein 
tensor 

V 
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(5-7) 



^ 6^ R . 

V V 2 V 

The Einstein field equations are then formed by equating 
the mixed Einstein tensor to another tensor kT^ called 

V V 

the mixed energy momentum tensor, an object vjhlch is related 

to T and is discussed later, 
yv 

Thus, the mixed form of the Einstein field equations 
as they are used in this V'jork are 

kT^^ = g\ (5-8) 

where 

^ . 



B. CONSTRUCTION OP THE CHRISTOFFEL SYMBOLS OF THE SECOND 

KIND FOR THE MODEL 

The task of constructing the field equations for the 
proposed model itself may nov/ begin as outlined in the 
preceeding section. From the field equations, the moment 
hierarchy will be determined. 

An important object that is immediately needed to form 
the Christoffel symbols is the second rank contravarlant 
metric tensor g^^. This object can be computed from Eqs. 
(i|-8) and (5-3). The nonvanishing elements of this tensor 



are 



gOO = 1 



S 



n 



-M 



g 



22 



= - e 

-M 

-e 
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(5-9) 

cont . 



g 



33 = 



-e 



-M 



cont . 
(5-9) 



sin^e 

The next step is to substitute the metric tensor of Eq. (^-8) 
into Eq. (5-1) which yields [yVj3]j the Christoffel symbol 
of the first kind. Because the second rank contravariant 

the 

Christoffel symbol of the second kind, by means of Eq. (5-2) 
is simplified. 

The calculation of is still somewhat lengthy but is 

entirely straightforward. The nonvanishing components of 

turn out to be 
yv 



metric tensor g^'^ is diagonal, the calculation of , 



ru 


= 


1 M 
f e“aoM 




= 


1 o M 

i r^e 3 oM 


0 

F33 


= 


^ r ^sln ^0e^3 


Toi 


= 


Fii = 1 


ru 


= 


^ 3 M 
2 r 


1 

r^ 


= 


fii = 1 9qM 


1 

T31 


= 


= 1 3 / 


1 

^22 


= 


- 1 + |) 


1 

T33 


= 


- r ^sin^e (3 M 

2 r 


2 

ro2 


= 


f20 “ ^ 9 0 



2 . 
r ■ 



(5-10) 
cont . 
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II 


1 1 cont. 

“ 2 T ^0^ (5-10) 


II 


r.I = 1 (3j.M + f) 


2 

T22 = 


1 


II 


T^ = i a .M 
^ $ 


= 


- 1 sin 20(2 cote + 9gM) 


fol = 


r* = 1 


Til = 


- I 3*M 


Til = 


fi3 = 1 ( + |) 


II 


^ sin^e ^ 


II 


fa = 1 (2 cote + 9 qM) 



r,! = I . 

C, CONSTRUCTION OF THE RICCI TENSOR FOR THE MODEL 

The next step is the calculation of the Ricci tensor 

R from Eqs. (5-^) and (5-10). Hov.’ever, it is anticipated 

p u 

that only the diagonal elements of the Ricci tensor are 
needed. (This v;ill be justified in the subsequent section.) 
These diagonal elements are 



Rot) “ 


1 a.^M + 1 O.M)' 


II 


+ i — i s'h oont^' 

2 ^2 e 2 ^23ij,2Q b cont. 

33 



2 






i + 1 3^M + J--H -77 

r r^sin‘^6 ^ 



R22 = - -^ r^e^a^M + i r^a^M + a^M 

2 0 2 r 6 



* i :T^ ^ ? r^(S^M)^ + 1 O^M)“ 



sin 6 



sln^e 



+ I cote a^M + I ra^M - ^ r2e*^(a(,M)^ 



R33 = - I r^sln^Oe^a^M + | r^sln^Sa^M + i sln^ea^M 

+ a|M - -|- r ^sin^6e^( a jM) + sln^S cot6 agM 



+ j r2sln20(a^M)^ + J- sin^eOgM)^ 



+ |- r sln^ea^M . 



(5-11) 



D. CONSTRUCTION HP THE DIAGONAL ELEMENTS OF THE MIXED 
EINSTEIN TENSOR FOR THE MODEL 

In order to construct the mixed Einstein tensor, the 
two objects R^ and R must be found. 

V 

The mixed Ricci tensor R^ is easily calculated from 

V 

Eq. (5-5)- Since is a diagonal tensor, R^^ can be 

written as 

R^ = (5-11) 

V ( X Jv 

/ 

where the parenthetical indices indicate no summation v;lth 
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respect to Equation (5-11) justifies the assertion that 

only the diagonal components of are needed if only the 
diagonal elements of are to be calculated because no 

off-diagonal mixing occurs by virtue of the dlagonality 
of g 

The Ricci scalar R is easily calculated from the diago- 
nal elements of R^^ by Eq. (5-6). Finally, the components 
of the mixed Einstein tensor given by Eq. (5-7) turn 

out to be 

G\ = - I (3oM)" + e”^[9"M + 9^M 

r 2 o 



+ ^ 

2 - 2 «d) ^2 D rr 

r^sin‘^0 ^ r 



cot 0 



2 '9“ ' If 2 • 2 o 

r r sin 0 



O^M)' + ^ (3^M)'] 



- 9oM - J OoM)' + I 3^M 

r 



+ i i 3 Jv l + - 3 M + y 3 M 

2 r^sin^0 ^ ^ ^ 2 ^2 0 



+ J Oj,M)'] 



^ (9oM)^ + e ^ 3 "m + ^ ^ 



2 r 



3?M 



2 r^sin^0 



r r 



9!m 



- I (9^M)“ + e-'''[ I 3jM + I ^ 3e” 



(5-12) 

cont . 
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I 



I 




1 



cont . 
( 5 - 12 ) 



+ i i 

2 r 



a M 



+ 



1 

? 



sin^e 



O^M) 

4> 



The off-diagonal components are not given because they 
are not needed. Hov;ever, calculations similar to the ones 
required to obtain Eq. (5-12) reveal that the off-diagonal 
components of are in general non zero and contain mixed 

V 

partial derivatives of second order and second degree prod- 
ucts of first partial derivatives v;ith respect to all four 
coordinates. 

It should be noted that if M is a special function of 
only the coordinates x° and r, namely 

M = G(x° ,r) (5-13) 

where G is given by Eq. (^-^), then the line element is 
identical to the line element of uniform cosmological m.odels 
[1]. In this case, as should be expected, Eq. (5-12) is in 
fact identical to the diagonal elements of G^^ for these 
cosmologies . 



E. CONSTRUCTION OP THE DYNAMIC EQUATIONS FOR THE MODEL VIA 
THE EINSTEIN TENSOR 

It should be noted that only tv;o "dynamic equations" 

for g are used to describe uniform cosmological models[l]. 
y V 

But each expression of Eq. (5-12) represents one side of 
four "dynamic equations." If the proposed model of this 
study is to pass to a uniform cosmological model for arbi- 
trarily small perturbations, then the number of equations 
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which should be treated as dynamical equations ought to be 
two rather than four. 

In order to see hov; to construct two equations from 
equation set (5-12), it should be observed that contains 

second partial derivatives with respect to all four coordi- 
nates x°,r,0, and (p . But each of the remaining objects 

contain second partial derivatives v;lth respect 
to only two of the three possible spatial coordinates. 

This assymetry can be removed if G^^^ is taken as one side of 
one of the desired two 'dynamical equations" and one third 

the sum G*^ . such that 
J 

G^ . = G* + G^, + G^ (5-1^) 

J 1 2 3 



is taken as the other. These two equations contain dynam- 
ical information concerning g and contain second partial 

yv 

derivatives with respect to all coordinates. They are also 
found to degenerate into the usual two cosmological equa- 
tions when G is replaced by M [1]. 

By virtue of Eq. (5-8), the two dynamic equations for 
the model are 

kT\ = G\ (5-15) 



and 




3 



(5-16) 



From Eq. (5-12), G^ . is given by 



= _ 332m - ^ (9 M)' + ^ 9^M (5-17) 

J 0 0 r ^2 e cont. 
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+ ^ O/i)' + O M)' 



sin^e ^ 



+ ^ O.M)^ + I (a^M) + 



sin^e ^ 



r r 



^ 3„tO. 

^2 y 



For convenience, G« can be expressed In terms of . 

0 J 



G»^ = + 39"M + I (a^M)^ , 



cont . 
(5-17) 



(5-18) 



F. CONSTRUCTION OF THE ENERGY-MOMENTUM TENSOR FOR THE MODEL 

To proceed further, it is necessary to express the 
energy-momentum tensor in a form that is specific to the 
metric used for this model. It is to be understood that 
the mixed energy-momentum tensor T^^ is also a stochastic 
process' because it depends, via the Einstein field equations, 
upon the metric tensor which is a stochastic process. The 
general form of the energy-momentum tensor is known to 
depend upon three quantities, D the proper mass density, 

P the proper pressure, and u the velocity of the medium. 

These quantities D and P are termed "proper" because they 
are referred to the proper frame which is that frame that is 
instantaneously moving with the medium. Prom any other 
frame, the well known effects upon the mass and volume due 
to the relative velocity u must be taken into account. 

The representation of the general energy-momentum tensor 
with respect to the metric that is used in this model is 
derived in Appendix D. All of the elements of this tensor 
are given by Eq. (D-19). The two related objects v;hich are 
important here are T®^ and T*^j ; 
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0 



TO 



= D + 




= - Du^ - — (3 - 2u^) (5-19) 

where 

D = proper density D x — - — 

1-u^ 



P = proper pressure P x (5-20) 

l-u' 

u = local velocity of the matter in units of c. 

Because of the recent discovery of a possible 3°K black- 
body background radiation [6, 23], a relevant energy-momentum 
tensor is one which describes disordered electromagnetic 
radiation of this type. However, this radiation is believed 
to be no longer Important to the dynamics of the universe 
(because now the matter density far exceeds the radiation 
density) although it is believed to have been in the more 
distant past. The energy-momentum tensor which is relevant 
to the current stage of expansion is one v.'hich describes the 
random motion of matter in the form of galaxies and clusters 
of galaxies. Eq. (D-M of Appendix D gives for these two 
cases 



P = ac^D 



( 5 - 21 ) 



where a = 0 if universe is matter dominated and 

a = 1/3 if universe is radiation dominated. 

Hence, by Eq. (5-20), 

P = ac^D. (5-22) 
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VI. MOMENT EQUATION HIERARCHY 



A. AVERAGE OP THE DYNAMIC EQUATIONS 

Equations (5-15) and (5-16) are the "equations of motion" 

for the assumed metric tensor of Eq. (4-8) and the energy- 

momentum tensor appropriate for that metric. These equations 

should be viewed as the analogue of the Navier-Stokes equa- 

tlon for V and P given by Eq. (A-1) because, as yet, the 

decomposed metric tensor of Eq. (4-9) has not been Inserted 

into them. -- 

Hov;ever , the substitution of the decomposed metric 

can be accomplished rather easily. A quick inspection of 

and G°^ reveals that they are made up exclusively of 

partial derivatives of the stochastic process M and a factor 
— M 

e . Prom Eqs. (4-10) and (4-12), it is easily found that 
e^ = e*^(l + t) (6-1) 

Q 

where e is associated v/lth the average metric and t with 

the fluctuating part of the metric. 

Equation (6-1) can be used to separate G*^ . and G°jj into 

J 

two parts, one v/hlch contains only the nonrandom function 

G and the other which contains the random function t in 

combinations of G and t. The part which is made up only of 

G terms must be G*^ . and G°g based on the average metric 

(because this part cannot depend on t) while the remaining 

part must be the source terms (because these terms must 

0 1 

depend on t). The resulting expressions for G ^ and G . 

J 



are (using Eq. (4-^) which is the definition of e^ and with 
R equal to and Eqs. (5-15) and (5-16)) 



= G\ 



2K 

R2 



3R 



R‘ 



3R' 

R 






1 + T 



3 illlL 

" U+t)" 



£f 

1 + T 



1+T 



3 

IT 2 

^ (1+t) 



1 + T 



3 



2 n 



(1+t) 



r ^sin^0 



2 -\ 



!ii 3 



cot 6 


f^e^l 


+ 2 1 


f \ 

^r^ 


r ^ 


1 + T 


^ . kr^ 


[i+tJ 



3k 

r2 



1 + T 



4r^ 

0 



( 6 - 2 ) 






. = G® 

J 0 



6r 



R 



39;t O^t) 



6R 

R 



V 

1 + T 



(6-3) 



If these equations are averaged, the result is the 
desired dynamic equations for the average variables <2^^^ 
and <T®„>. 

These averaged equations can be Immediately simplified 

by using the assumption of spatial stationarity that is 

discussed in Ch. IV. Thus, all averages of the form 

9.<h(T)> vanish by virtue of Eq. (4-20). Therefore, in 
0 



particular , 



<9 . 
0 



1 + T 



> = 0 



Prom this follow the identities 

2 



(l+T) 'J 



> = < 



9 !t 

.1 



(1+t) . 



and 



V 



L(1+t) 



> = 0 



Thus, with Eqs. (6-5) and (6-6), the average of Eqs 
and (6-3) are 

k<T°„> + 



3R’ ^ 




> + ^ < 


■(9oT)'' 


R ^ 


1 + T^ 


> i- ^ < 


.(1+t)'. 







r 2 1 




r 

(90t) 


1- e-®- 




(9rT) 


> + < 


Tf ^ 


s, 


.(l+T)^ 


r=^ 


.(1+t)'.^ 



1 






. _ 3k ^ 


f N 

X 


r^ sin^ 6 


.(1+T)'. 




R2 


1+T 



3k 3R' 



R‘ 



R' 



5 <tJ > + ^ <[!ii 

O ^ ^ R 1+T 



[3.^- 


> - ^ < 


’(9ot)'' 


Li+tJ 


> - IT < 


L(1+t)'. 



+ < 



■(92t)' 


> e-®. 




r 2 1 


.(1+t). 


12 


V. 


.(1+t)'. 



+ 1 - < 
r^ 



(1+t) . 



> + 



r ^sin^O 



( 9 , T ) 
.(1+t) 



2 1 



(6-H) 

(6-5) 

( 6 - 6 ) 

(6-2) 



(6-7) 



( 6 - 8 ) 
cont . 
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are 

Eq. 



ll+-r) R^ R 



( 6 - 8 ) 



If Eqs. (5-19) which give expressions for t\ and T j 
used, Eqs, (6-7) and (6-8) become (with the aid of 

( 5 - 22 )) 

3R * 1 ° 

k<D( 1 + ctuM> + "X" 




r ^ sin^ 0 



_ 3k _ 3R1 



R^ R^ 



3R' J ° 

J< <Du'(2a-l) - 3aD> + — 

3 






.(l+t) 



(6-9) 
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Ogt) 
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(3 t) 

■ q 
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, 3 

,(1+t) - 


> ^ 

r^sin^6 


1(1+t) - 


1 
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R^ 



2R” _ R_L 
R R^ 



(6-10) 



,tI 0 ns are the desired dynamic equations, 
,e two equations 
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B. CONSTRUCTION OP THE SPECIAL MOMENT EQUATIONS 



As stated in Ch. II , it is necessary to determine the 
functional form of the extra source terms of Eqs. (6-9) 
and (6-10) in order to solve for the average quantities R, 
<T° „>, and <T*^ .>. 

The basic idea is to construct additional equations 
for the unknown correlations. This is done by first multi- 
plying the unaveraged Eqs. (6-2) and (6-3) or combinations 
of these by some "special functions" which involve pov.'ers 
or derivatives of x and then averaging them. As might be 
expected, the nonlinearities of the unaveraged equations 
cause more unknovm correlations to be formed v;hich requires 
the formation of more equations. Although this process can 
be carried out ad Infinitum, a special truncation of this 
process is required which imposes an additional restriction 
upon the model. The truncation that will be used is the 
higher-order-moment discard approximation that is discussed 
in Ch. Ill in connection with the fluid turbulence hierarchy 
This means that the higher order moment that results from 
the nonlinearity is discarded because of smallness. Hence, 
only the first "generation" of the moment hierarchy is 
needed and it turns out that three special equations are 
sufficient (see Ch. Vil). 

The two special combinations of the unaveraged dynamic 
Eqs. (6-2) and (6-3) which are used as the basis for develop 
ing these three equations are 
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2 



( 6 - 11 ) 



and 



k(T®o + T-^ .) = g\ + . 

J J 



= 2G 



6Rt. 33 qT _ 2 



R 



1+T 2 



(1 + T ) 



6R 
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1+T 



^(TO IJ.J ) = G® - G-J . 



6rm 33oT 3 6R' 

R 1+T 2 , s 2 R 



^0^ 



1 + T 



(6-12) 



(1 + T ) 

If Eq. (6-11) is multiplied by the "special function" 

1+T and averaged, the first moment equation is formed. The 
resulting equation can be simplified by again using Eq. 
(4-20) which implies that 

3?<h(x)> = 0 . 

J 

If h(T) is taken to be ln(l+T), then the following identity 
results 

(s.t) 

(6-13) 
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.1 + T . 


> “ 


.(1+t)\ 



Thus, the simplified first moment equation is 
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If Eq. (6-12) is multiplied by t and then 9 qt, the 
second and third moment equations are produced. 



(c<(T° - Jt> - 3< 
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■(3^t) (t)’ 
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(6-15) 
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Again, using equation set (5-19), these three equations 
become 

K<D(l-u2)(i-3a)(l + T)> - I 



+ 1 - < 



- ( 1+t) 
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K<[D(l+3a) + Du2(l-a)]T> - 3< 
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(6-19) 
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VII. TRUNCATION OF THE MOMENT EQUATION HIERARCHY 



In both Eqs . (6-9) and (6-10) it can be seen that the 
source terras contain (1+x)^ in their denominators. Because 
(1 +t)'^ can be expanded in a power series in t, the source 
terms can be expressed as an infinite sum of moments of 
various orders in x and its derivatives. As mentioned 
before, an equation for each of these moments can be gener- 
ated by multiplying the unaveraged Eqs. (6-2) and (6-3) 
by an appropriate quantity and averaging. The nonlinearities 
of these equations create a dependence of the desired moment 
on higher order moments. Because this occurs for all moments, 
a method is needed to truncate the system in such a v;ay 
that only a finite number of moments and required equations 
in the hierarchy are left. In Ch. Ill three possible trun- 
cation schemes are presented. In this v/ork, the higher- 
order-moment discard method is chosen by assuming that 
most of the probability mass of the random variable x(x^,0 
is concentrated around x = 0, so that the higher order 
moments can be neglected. This obviously limits the appli- 
cability of the theory to models of small nonuniformity 
for which |x|<<l (see Ch. IV-B) . 

However, care must be exercised in affecting this trun- 
cation. It is not correct to just replace (1+x)^ by one in 
each of the denominators of the above mentioned equations 
just because this is the first term of the expansion. This 






can be seen by an example. 

f \ 



< 



Ii+tJ 



> 



This term can be 



Consider the typical term 
expanded as follov;s 



< 




O^T> - <(8|jT) (t)> + . . . 



(7-1) 



Since O x> = 3 <t> = 0. 

0 0 ^ 



< 




< ( 3 jT ) ( T ) > . 



(7-2) 



However, if (1+x) were falsely replaced by one, the follov;- 
ing incorrect conclusion would result which disagrees v;ith 
the correct Eq. (7-2) 



< 




<3 „x> 
0 



(7-3) 



Hence, an expansion must be made for every term and that 
lowest order moment term v/hich does not vanish by virtue 
of having' a zero average is assumed to be a "good" approxi- 
mation to the term itself. 

With this "weak fluctuation" assumption, Eqs. (6-9) 
(6-10), (6-17), (6-l8), and (6-19) become 



K<D(l+au^ ) > 



3 ^ 
2 R 



3 <x^> 
0 




2 

t) > 





3rT 



2 -I 2 

> + — <(3^x) >+ 

■K. 2 



0 



> 

r^sin^e 



+ 3k _ 3k _ 3^ 

R^ R^ • 

^ <Du2(2a-l) - 3aD> - | " I 

3 ^ K « ^ cont. 
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+ f <Oo'^) ^ ~ ll ^ ‘^|<(9^t)'>+ ~ <Oe'^)'> 



1 2 

+ — > 

r^sin^e ^ 



■ + — <T ^ > 

R^ 



k 2R" R 



R^ R R- 






r^sln^S 



cont . 
(7-5) 



<O^T) > 



+ I <<D(l-u')(l-3a)(l+T) 

_-15k 15R'" 15R” 

2 “ ; (7-d) 

R^ R^ R 

K<[D(l+3a) + Du^(l-a)]T>- |- 3p<T^> 

QR I , 2 

- > = - 3<Oo'^) > (7-7) 



^ IJP t 2 

9o<(3o'^) > + ^ <00-^) > 

= I K<[D(l+3a) + Du^(l-a)]3jT> _ (7-8) 

In the limit of zero fluctuations, Eqs . (7-^0 and (7-5) 

become the usual tv;o cosmological equations for p and R 

that are given by Adler [1]. Thus, it is natural to look 

upon these equations as dynamic equations for p and R. 

The source terms v/hich appear in these equations contain 

2 2 

the three functions <x^> , <(3 (,t) >, and (<(9 ^t) > 

+ ^ + F2 ~ sfn ~ 0 
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equations are needed to determine these variables v;hlch 
justifies the statement made in Part B of the proceeding 
chapter that only three equations of the hierarchy need 
to be calculated. These three equations are Eqs . (7-6), 

(7-7), and (7-8). Although the above functions <t^>, 

2 2 1 2 1 2 
and {<(9 t) > + — ^<(9 t) > + - <(9,t) >} 

appear in these equations, additional moments also occur 

which relate the density D to t and D to 9pT. It is shovm 

in the succeeding chapters that these moments can be related 

to the existing moment <x^>. Thus, Eqs. (7-^) to (7-8) 

are soluble among themselves. Because Eqs. (7-6) to (7-8) 

govern the dependence of the source moments on R, they 

are called grov^^th equations to convey the idea that they 

determine the development or growth of a perturbation t v/ith 

increasing R. 
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VIII. REVIEVJ OF THE RESTRICTIVE ASSUMPTIONS 



Before proceeding further, it might be helpful to reviev; 
the restrictive assumptions that are used to obtain Eqs. 

(7-^) through (7-8). The first is the very restrictive 
form of the metric tensor. This form excludes many types of 
gravitational wave phenomena. All that is seen by an observ- 
er would be a local "expansion" or "contraction" of space 
because, as stated in Ch. IV, it would alv;ays be possible 
to construct an expanding or contracting Robert son-V.'alker 
geometry that is "tangent" to the given space at any given 
point in space-time. The exclusion of such v;aves must 
certainly exclude any hope of using this model to describe, 
in a realistic way, the violent, multiply connected space 
that Wheeler envisions must typify the quantum fluctuations 
of space-time. 

The second and third restrictions are the restrictions 
concerning the Robertson-Walker form of the average metric 
and the spatial statlonarity of the fluctuating part. These 
two assumptions rely more on intuition than upon scientific 
rigor. They imply that there exists some sort of spatial 
uniformity in the statistics of the assumed fluctuating 
geometry. Hov/ever, it is to be understood that a "new" 
cosmological principle is not being seriously proposed here 
by the use of these two assumptions. That these assumptions 
are used at all reflects a need to simplify the mathematics 
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of the moment equations in some reasonable way by appealing 
to some plausible symmetry arguments. Thus, these tv/o re- 
strictions should be viewed only as steps in this particular 
problem whose purpose is to demonstrate a technique. 

Finally, the fourth restriction is the assumption 
that the higher order moments are negligible ov;ing to an 
assumed smallness in the anisotropy and inhomogeneity. This 
assumption is a restriction upon the applicability of the 
model. That is, in those solution domains for v;hich the 
magnitude of t is comparable or larger than one, the use of 
the model is unwarrantable. If large amplitude fluctuations 
are to be considered, then perhaps a truncation based upon 
a phenomenological or a cumulant-dlscard approximation would 
be useful. However, in the case of the latter, it might 
be assumed that the probability distribution of the random 
variable t is normal. Then, the higher order moments can 
be shown to be products of the lov;er order ones [l8]. Such 
an assumption is made in fluid turbulence theory [3]- But, 
in the case of this model, x must be greater than or equal 
to minus one as given by Eq. (4-15). This creates difficulty 
with a normal distribution where the probability of x being 
+T is equal to the probability of it being -T v/here T is an 
arbitrary positive number. 



» 
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IX. SPECIAL UNIVERSE MODELS 



From this point on it is convenient to consider tv;o 
types of models which correspond to tv;o values of a . The 
first, for which a is equal to one third, corresponds to a 
universe made out of electromagnetic radiation while the 
second corresponds to a univers’e v;hich consists of cold 

p 

matter for which a is zero. (a is equal to . See Eq. 



It is not practical to treat these models simultaneously 
by leaving a arbitrary. The cold dust model can be solved 
exactly owing to the fact that all matter moves freely along 
geodesics without being accelerated by any pressure gradi- 
ents (pressure is zero everyv;here ) . But the radiation model 
does involve pressure gradients which do not allow the 
elements of the radiation fluid to move freely in space [9] 
Therefore, these two types of models are sufficiently differ 
ent from one another to merit separate consideration. 

A. RADIATION FILLED UNIVERSE MODEL 



For the radiation filled universe, Eqs. (7-^), (7-5), 
(7-6), (7-7), and (7-8) become 



( 5 - 22 ).) 





cont . 



(9-1) 



5 *^ 



+ » = - a . ai 

R2 



cont . 
(9-1) 



_ -'le. _ 3 R_L 



2 R °o 



f + F <(3.t) 






e"^ 
12 ^ 



<(9j^t) > + — <(9qT) > 

2 V 



<(9 , t ) > [■ + — < T ^> 



r^sln^e J r2 

_ k 2^' _ RJ_^ 

R2 



R R2 

^ e *^ i <(8 t ) > + — <( 9 ^ t ) > + — 



(9-2) 



r ^sln^e 



< ( 9 , T ) > 



-15k _ 15R' _ 15R" 

R^ R2 R 



(9-3) 



2K<Td> - I 9 ^<t^> - 

= - 3<(9ot)"> (9-4) 

2 ilR t 2 U 

where 

p = <D(1 + ^)> (9-6) 

and 

d = D(1 + 1^) - p (9-7) 

It should be noted that P Is equal to the average of 1°^^. 
Adler [1] points out that 1°^ can be interpreted as the 
total mass density with respect to the given coordinate 
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system and therefore includes relativistic corrections to 
the proper mass density. 

1 . Elimination of the Spatial Derivatives 

Now consider the terms v/hich occur in Eqs. (9-1) 

2 

and (9-2) of the form <(3 ^t) >. It should be noted that 
only and 9^ are partial derivatives with respect to a 
length (in the time and r directions). The remaining par- 
tial derivatives 9g and 9^ are partial derivatives v/ith 
respect to the coordinates 0 and 4> which are not lengths 
in themselves. An inspection of the line element of Eq . 
(^-2) suggests that 9^, e~*^^^9^, e~^^^ 1/r 9g, and 
l/rsln6 9^, all of which have dimensions of inverse lengths, 
play the same roles in the x“,r,0, and (fi directions respec- 
tively. Because disturbances of the metric propagate with 
the speed of light [11, 32 ], it v;ould be expected that the 
characteristic Inverse lengths associated with these quanti- 
ties are comparable in magnitude. Thus, it seems reasonable 
to expect that the spatial terms 



r 2 1 n 2 -G 

e"‘'<(9 T) >, e ^<(3„T) >, — £ 

^ r^ (r sin0) 






could be expressed as mere repetitions of the temporal term 
2 

<(9^t) > multiplied by a coefficient of order one. 

Such a fact would simplify Eqs. (9-1) and (9-2) 
because the terms involving the spatial partial derivatives 
would be changed into the already existing terms Involving 
the time derivative. 
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This consolidation of the terms is performed in 
Appendix F. Equation (9-3) is used to eliminate the spatial 
partial derivative terms from Eqs. (9-1) and (9-2). An 
appropriate renormalization is then performed to bring the 
right hand sides of the resulting equations back to "canoni- 
cal form" as they appear in the right hand sides of Eqs. 
(9-1) and (9-2). Equation (9-^) is used to change the form 
of the "renormalized" equations into a more convenient form. 
This last step is not part of the consolidation itself 
since it just allows the resulting equations to be presented 
in a more desirable form. The transformed pair of Eqs. 

(9-1) and (9-2) which are void of spatial partial deriva- 
tives are 



Kp 



3 R' 



Ij- K<Td> + i;- <(3 „t) > + B" r" 



3 IL. 

^ R2 



<T^> = 



3k _ 3^ 

R2 R2 



(9-8) 






1 |1 



R^ 



k 2R" R 



(9-9) 



R^ R R- 



2 . Inclusion of the Dynamics for the Radiation Field 
Because Eq. (9-3) has been used to remove the un- 
known 



-G 



2 1 .2 1 , . 2 
<(3 t) > + — <(3^t) > + <(3,t) > 

r j,2 e - 



r^sin^S 



from Eqs. (9-1) and (9-2), there exists now only four 
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equations In the system namely, Eqs. (9-^), (9-5), (9-8), 

and (9-9). But this system of equations is insoluble because 

they represent four coupled nonlinear differential equations 

in six unknovm functions of the time coordinate x° . These 

2 

unknown functions are p, R, <t^>, <(9gT^) >> <Td>, and <O^T)d>. 

However, there is a way to express <xd> and' <(9gT)d> 
in terms of the other four unknovm functions. It is not 
unreasonable to expect that there should exist some rela- 
tionship between the density fluctuation d(x^) and geometri- 
cal fluctuation t(x^) because in some sense the geometry 
fluctuations "drive" the density fluctuations. But it must 
be remembered that the radiation fluid possesses internal 
pressure forces which push the fluid elements about as 
well as do work against any expansions or contractions of 
space which might occur v/lth amplitude t [9]. Thus, the 
problem of determining this relationship betv/een d and t is 
in general very complicated. The correct v.'ay to do this 
is to use the equations of motion for the radiation fluid 
which are given in tensor form by [1] 

T^^; V = 0 (9-10) 

where 

T^'^; V = 3 T^*^ + f^ T*^'^ . 

’ V va av 

Indeed, the relationship betv.'een d and t is approximated 
in Appendix I via these equations and is expressed in terms 
of the correlations <xd> and <(3jjT)d>. 
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An instructive yet simpler v;ay to determine the rela- 
tionship between d and t is the following. Consider a small 
3-volume element V at rest with respect to the coordinates 
(r, 0j (}) ) . Let M be the total relativistic mass content 
of this 3-volume. The total relativistic mass density D 
is thus given by 




A small change AD of D can be in general a result of a small 
change AM of M and a small change AV of V. That is. 



AD 



m 

V 



MAV 

V2 



(9-12) 



Now the quantity d‘*x is known to be a scalar 

invariant v/here g is the determinent of the matrix g^^ 

[1, 32 ]. In particular, from special relativity, this quan- 
tity is the product of the proper time interval and the 
proper volume. Because x® has already been identified as 
the time coordinate (see Ch. IV), the volume V must be 



V = d®x . (9-13) 

In terms of the coordinates r, 0, <p , 

V = drdOd(}) . (9-1^) 

The use of Eq. (4-8) gives 

V = e^^'"^'^^(dr ) (rde) (r slnedq). (9-15) 

But e^ = e*^(l+'^) (9-l6) 
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so that V = V(l+x) 



(9-17) 



3/2 



where 




(dr)(rd0)(r sin6d({)) . 



The change of mass AM can occur because of the work 
done by the pressure against the change of volume and by 
the flux of mass out of the volume. If this flux is tempo- 
rarily ignored (i.e. all mass moves along geodesics), then 
AM can be expressed as 



where P is the total relativistic pressure. The relation 



AM 




(9-18) 



P=io^D 



(9-19) 



Implies that 



AM 




( 9 - 20 ) 



Then Eq. (9-12) becomes, using Eq. (9-20) 



AD 




(9-21) 



But Eq. (9-17) with small becomes 



V = V + I Vt . 



( 9 - 22 ) 



If AV is replaced by ^ Vx , V by V, AD by d, and D by P in 
Eq. (9-21), then this equation becomes 



d 



2px 



(9-23) 



Hence, the correlation function <xd> and <(9pX)d> become 



<rd> = - 2p<x^> 



and 



<Oj,T)d> = - p9q<t2> 



( 9 - 2 ^) 



In Appendix I, the effect of "mass leakage" from the volume 
is approximated by means of tv;o variable parameters X and 3. 
It is shown that a reasonable form for the correlations 
<xd> and <(9^T)d> are 



The coefficients X and 3 may be thought of as representing 
a zeroth approximation of those functions of x°, which 
correct the zero leakage correlation Eqs. (9-2^). These 
coefficients, then, become open parameters which character- 
ize a given model. Certainly, if mass leakage were not 
important, then X and 3 would both be one. 



Eqs. (9-8) and (9-9) and into the growth Eqs. (9-^0 and 
(9-5) gives 



<xd> = - 2Xp<T^> 



<(9pT)d> = - 3p9o<T^> . 



(9-25) 



The substitution of Eqs. (9-25) into the dynamic 




2 




3k 3R' 



( 9 - 26 ) 




2 




k 2R” R' 



(9-27) 



r2 R R^ 



6l 



aJ<T^> + + I Xkp<t^> = 2<(3^t)% (9-28) 



2 hXD f 2 

<(3(jT) > = 



BKp9n<T^> 



(9-29) 



3 • Solution of the Moment Equations which Govern Grov;th 
of the Fluctuations 

The four Eqs. (9-26) through (9-29) now contain the 

2 

four variables p, R, <t^>, and <(3 (,t) > and the system Is 
now closed. However, the solution of these four equations 
Is not easy because they are nonlinear and nonlinear prob- 
lems do not In general yield easily to analysis. But It 
Is possible to take advantage of the requirement that the 
moments v;hlch appear In Eqs. (9-26) and (9-27) be small 
compared to icp. That Is, the additional stress terms cause 
a small perturbation of the usual, fluctuation free cosmol- 
ogical solutions for p and R which are well known. These 
solutions, derived In Appendix E, are 

3a„ 



Kp = - 



(9-30) 



R‘ 



and 



2 a 

R' = - k 

R^ 



(9-31) 



where a Is a constant of Integration and k Is the curvature 
r 

characteristic and may take one of three values -1, 0, +1 
(see Ch. IV). 

Nonlinear problems can In many cases be simplified 
by the replacement of a variable with a very close approxi- 
mation to that variable [30]- such a replacement Is 
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made for R' and p in Eqs. (9-28) and (9-29) with the fluc- 
tuation free solutions Eqs. (9-30) and (9-31) and if all 
derivatives with respect to x° are transformed into deriva- 
tives with respect to R by the chain rule for differentia- 
tion, these two equations become soluble between themselves. 

However, as can be seen from Eq. (9-31), the effect 
2 

of k upon R’ is not significant unless R is large enough 
so that a^/R^ is of the same order of magnitude as k. It 
is well understood that that value of R for which the effect 
of k is important to the dynamics is the current value of 
R. For this reason, much research in observational cosmol- 
ogy is devoted to the determination of k on the basis of 
data on current expansion characteristics [25]- But the 
applicability of the radiation model is only valid up to 
the time when the radiation density is equal to the matter 
density for which the associated R-value is much less than 
the current value of R by several orders of magnitude [8]. 
Thus, if the R-domaln of the solutions <x^>(R) and 
<(8qT)^>(R) Is restricted to values of R sufficiently less 
than the current value of R so that a^/R^» |kj , then the 
following substitutions for R' and p in Eqs. (9-28) and 
( 9 - 29 ) are valid. 




R“ 



and 




(9-33) 
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It should be noted that the neglect of k in cosmol- 
ogical problems is not nevj. k is neglected, for example, 
by Hoyle and others in their expression for R ' for exactly 
the same reason that it is neglected here [10]. 

The growth equations become, using the above approx- 
imations 

R23^<t"> + R3j^<t"> - 8x<t"> = (9-3'^) 



’R 



. 2 ' 
< (3oT) > 


+ i 


' 2 ^ 
<(3oT) > 


= M 


^r 

^ J 


^ R 


^r 


R“ 



’R 



( 9 - 35 ) 



These equations are solved in Appendix J for <t^> 



and <(3 |jT) >. The solutions are 



<T^> - F^/3 + R + R 



( 9 - 36 ) 



< (3 qt) > _4 i -j 

= - “iXP./jR + I1/3R + 



( 9 - 37 ) 



where i = - ^ + 2/2 ( A+B ) 

0 = - - 2/2(A + B) 

and > ^1/3 ’ ^1/3 arbitrary constants of integration. 

1] . Special Submodels 
a. Discussion 

In general, a given combination of the parameters 
I ^^3 , , and ® noriuniform model which 

is difficult to treat analytically. Yet, there are special 
values of these parameters which designate models that are 
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relatively easy to treat and which are interesting. These 
are the three models for v:hlch only one of the three arbi- 
trary parameters are non zero. That is, 



< ( 3 0 T ) > . . 

(1) <T = > = P^/3 ; = - 4XFJ/3R-'' 



( 9 - 38 ) 



^1/0 T-l-il ^(3nT)> J 

(2) <T^> = — r--- ^ • — ■“ T P 



(9-39) 



^1/0 i+il J 

i3) <T > R , — 



(9-^0) 



Actually, not much generality is lost by considering the 
above three models separately. The I 2/3 » ^1/3 ’ ^1/3 

terms of the general Eqs. (9-36) and (9-37) are functions 
of R to a power. In general, liov/ever, these powders are 
quite different from one another (i.e. try A, B equal to 
one which are the zero leakage values for A and B). Because 
R changes by many orders of magnitudes during expansions 
and contractions, only one of the terms v;ould dominate 
anyway (except for those "particular" values of R for which 
two or even three of the terms would be of the same order). 
Hence, most of the models with a given combination of I 2/3 ' 
^ 1/3 would degenerate into one of the 

above three models for almost all R-values. > 

However, there is a difficulty with the F, ,^- 



model. Equation (9-25) says 
<xd> = - 2Ap<x^> . 
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X must therefore be a parameter v/hich measures the response 

of the radiation fluid to a given geometric perturbation. 

It seems unreasonable to suppose that the fluctuation d of 

the total mass density is positive for positive t- That is, 

it seems that only models that are physically reasonable 

are the ones for which X is positive. But, from Eq . (9-38), 

only the negative X values permit the obviously positive 

2 

functions <( 8 jjT) > and <t^> to be positive (a^ is positive 
by virtue of Eq. (9-32)). Hence, the F^^^-model does not 
seem to have much physical meaning and is therefore not 
considered in this work. (If some positive X values v;ould 
be acceptable then there would be grounds for exploring 
the content of the P^^^-model.) 

The two remaining models, namely the 
D^^^-models, are very different from one another. The 
Il/ 3 -model represents a model which becomes more nonuniform 
with increasing R while the depicts a model for 

which the nonuniformity decreases v;ith increasing R. 

Because there is no intrinsic reason why positive X values 
are not acceptable, both of these models can be considered. 
(I ^^2 is used to denote the arbitrary constant of an 
increasing type perturbation field of the metric for a 
radiation, a equal 1/3 universe and used to denote 

the decreasing type. The same nomenclature will be used 
to specify the Increasing and decreasing type perturbations 
of the dust, a equal zero universe but with a zero subscript.) 
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It is of Interest to see hov; the fluctuation 



stresses which correspond to both the I 2./3 D^^^-type 

perturbations act upon the average total mass density p 
and the average scale length R. This can be seen if the 
stress terms which appear with the average density p in 
Eqs. (9-26) and (9-27) are expressed in a form which cor- 
responds to the functional form of the D^^^-type 

fluctuations given by Eqs. (9-39) and (9-^0). Equations 
( 9 - 26 ) and ( 9 - 27 ) become (with ^ <x > neglected to 
Kp<T^> for reasons already stated) 



( I -flue t nations ) 

KP . c(i^/3)r^ [60X . 6/m mi 

= _ 3k _ 3R_!_ (9-^1) 

R^ 

Kp xpi [12X + 8B - 2/2('A+'B)] 

"3 72 

= _ _ Rll (9_il2) 

R^ R R^ 



(D^y^-f luct nation s ) 



Kp + C(D^y^)R 



j [ 6 OX + 2^3 - 6/2(A+g)] 
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3k 3R' 

r 2 R^ 



- I p - 0 ( 0 ^/^)^' 



[12X + 8B + 2/2(A+g)] 

72 



2 

k 2^' _ L!_ 

R^ R R^ 



(9-^3) 



(9-^4) 
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9 I-| 3- 9 D-, 3, 

where and CCD^^^) - ^i-X. 



1 = - + 2/2(A + 6) and j = - 4 - 2/2(A+g) . (9-45) 



It should be pointed out that of the three 
parameters C(A), A, and 3 of Eq. (9-^5) (v;here A represents 
If/S ^1/3^’ Calends itself to physical measurement. 

It must not be forgotten that is proportional to a para- 
meter vrhich Is the result of an integration and therefore 
is an arbitrary constant vjhlch must be fixed by some obser- 
vation. But A and 3 are not the result of any integration. 
They are variable parameters which, in a sense, measure the 
mathematical uncertainty between tv;o pairs of correlation 
functions (see Eqs. (9-25)). If it were known hov/ to 
exactly determine these relations, then 3 and A v/ould be 
known numbers and hence be of nonvariable status. As it 
is, they have to exist as variable parameters. (In the 
dust model considered later, the dynamics of the matter 
are determined exactly. Hence, the correlation functions 

can be found precisely and A and 3 do in fact, for that 
case, become knovm numbers.) 

The solution of the dynam.ic Eqs. (9-^1) to 
(9-44) are given in the next section, 
b. Classification 

The question naturally arises: What is the 
effect of a given combination of fluctuation parameters 
C(A), A, and 3 upon p and R? It has already been pointed 
out in Section 3 of this chapter that the effect of the 
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fluctuation stresses upon the average density p and scale 
length R is small. Thus, the actual solutions are small 
perturbations of the fluctuation free solutions and it 
therefore seems appropriate to compare the solutions of 
Eqs. (9-^1) through (9-^^) for p and R with these known 
fluctuation free solutions. 

This can be done in the follo’wing manner: There 
are only four distinct vmys that the fluctuation stresses 
can perturb the solutions for p and R from their fluctuation 
free solutions. That is, at a given time the density 

p and the scale length R can each be perturbed above or 
belov; their fluctuation free values. Let each of these 
four type of perturbations be designated by the letters K, 

L, M, and N in the manner shown in Fig. 1 for the increasing 

and Fig. 2 for the decreasing D^^^-fluctu- 
ations. In both cases, a given letter indicates v/hether 
the p and R are greater or smaller than the fluctuation 
free solutions for p and R respectively. For example, a 
model for which p(x°) and R(x°) are both greater than their 
respective unperturbed values p*(x°) and R* would be clas- 
sified as type M. 

The manner in which the solutions p and R are 
perturbed from their fluctuation free values which are 
specified by the above simple classification scheme is 
determined solely by the two parameters X and 3* The mag- 
nitude by which a perturbation of a given type occurs is 
determined by the parameter C(A). Thus, X and 3 determine 
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the qualitative behavior of the model v;hile C(A) determines 

its quantitative behavior. 

In order to determine the classification of a 

model which is characterized by the three parameters C(A), 

X, and 6 where A refers to either the increasing I^^ 2 ~type 

perturbation or the decreasing D^^^-type, the solutions for 

p and R are needed. In Appendix L the solutions for tcp 
2 . 

and R* are given in terms of some special parameters f ^ , 
b^, and m^ v;hlch are functions of only the three parameters 
C(A), X, and 3. These solutions are 



-4 % 

Kp = - 3a^R - 



2 -2 
R’ = a^R ^ + b^R 



( 9 -^ 6 ) 

(9-^7) 



The unperturbed solutions p* and R* are those 

for which f^ and b^ are both equal to zero. It is easily 

seen that the direction of displacement of p from p* is 

determined solely by the sign of f^. The sign of the para- 

2 

meter b^, on the other hand, tells v/hether R' is larger 

2 » 

or smaller than its unperturbed value (R') for all values 

of R. If, in the case of the I^^^-type model, a forward 

time integration of Eq. (9-^7) for various values of bj 

is made from R equal zero, a family of solutions R(x°) 

can be generated. The displacement of any one of these 

2 2 ^ 

solutions from R‘^ depends on the ratio R' /(Jl ' ) • For example, 
if this ratio is greater than one for all R-values, then R 
is necessarily perturbed above R* because all solutions of 



1/3 
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the family are "pinned" to the same Initial value. Because 



■ 1 / 3 ’ 



the size of this ratio compared to one is determined by b. 

it follows that b-p also determines the direction of the 

^ 1/3 

displacement of R from R* . It can be shovm that the sign 
of bp^y^ 3-lso determines the displacement of the decreasing 
solution of R from R* if the Integration is begun 
at a sufficiently large R and performed backv;ards in time. 

In Appendix L, the signs of f^ and b^ (and hence the classi- 
fication letters) are determined for all values of \ and g 
for both the ^1/3 fluctuations. The results 

are graphically displayed in Figs. 3 and 

It might be asked how a given combination of X 
and 3 affect the exponent of R in the fluctuation stress 
term. The dependence of this exponent is presented in Fig. 

5 as a plot of the lines of constant j and i in X-B space. 

This concludes the basic discussion and calcu- 
lations for the radiation model. In the next section, a 
similar treatment is made for the dust model. The discus- 
sion and comparison of these two models are given in 
Chapter X. 



B. DUST PILLED UNIVERSE MODEL 

Some of the calculations and discussions in this section 
are similar to ones of the previous section which deals 
with the radiation filled universe for which o. is equal to 
one third. Here the model for which ot is equal to zero is 
treated. It is assumed that the reader is already familiar 
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with the previous section so that otherwise redundant expla- 
nations may be omitted. 

For the pressure free universe, Eqs. (7-^), (7-5), (7-6), 
(7-7), and (7-8) become 



OR! O 2 

Kp - 2 R“ TT > 



- ^ e ^•|<(9 ^t) > ^ — <(8,t) > 



r^sin^S 



+ 3k<T^> = _ 3k _ 3R1 
R2 r2 



(9-^8) 



- I|- - I I <0,t)'> 



^ e-° 

12 



<(a t) > + ^ 2 '!"7” 

^ r^ ® r^sin^e ^ 



^ k 2 k 2R" R' 

R2 



R2 R R2 



(9-^9) 



5 ^-G, 
TT ® ^ 



<0 t) 


\ . ^ 


X 




5 

J KP - 


K<Td> 


3 r\ 2 ^ 


2^ 3R' 


2 


^ R 



r^ sin^e 



15k 15R' 15R 



< (SATf ) > 

4) 



(9-50) 



R^ 



3 <T^> = 



R^ 



R 



- 3<0^t) > (9-51) 



2 ^2 2 
0 '’O’' R ' 0 






(9-52) 



where 



p = <D> 



and 



d = D - p 



(9-53) 
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It should be noted that for a equal to zero, the pressure 
is Identically zero. Hence, there are no pressure gradients 
available to accelerate the matter. This implies that 
matter, when left alone, v/ill follow geodesic lines. If 
these geodesics are taken to be geodesics of zero velocity 
then the matter comoves with the coordinate system (r,6,(})) 
and the above equations follov;. From Eqs,(9-53) it is seen 
that p is the average of the proper density D because the 
relative velocity u is zero in which case D is equal to D 
(see Eq. (5-20)). 



1 . Elimination of the Spatial Derivatives 

In Eqs. (9-^8) and (9-^9) j the terms v/hich contain 
spatial partial derivatives may be eliminated explicitly in 
a manner analogous to the v/ay they are removed from Eqs. 
(9-1) and (9-2). This computation is performed in Appendix 
G. The resulting two equations are 

Kp + I <(9oT) > + I f 

2 



3 k_ _ 3k _ 

^ R2 R2 R2 






2 

t) > 






12 



K<Td> 



1 k ^ 2^ _ k 2 ^ ^ 

-I- <T > - - - - 

^ R2 R2 R R2 



( 9 - 5 ^) 



(9-55) 



2 . Inclusion of the Dynamics for the Matter 

In order to solve the system of equations, it is 
necessary to express <xd> in terms of existing correlation 
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functions. This is also done in the previous section but 
is complicated by the existence of nonzero velocities. 

For the pressure free model this calculation is much easier. 

Consider Eq . (9-12). Because the velocities and 
the pressure are zero, there is no mechanism to change the 
mass content of an expanding volume. Hence, AM must be 
zero and Eq. (9-12) becomes 



AD = - 



MAV 



(9-56) 



With Eq . (9-22) for small Eq. (9-56) becomes 
(also with the usual definitions for the increments as in 
the previous section) 

d = - I pT . (9-57) 

Hence, <xd> and <(9oT)d> must be 

<Td> = - P<T^> (9_58) 

and 

<(3ox)d> = "I • (9-59) 

Equation (9-57) is also derived in Appendix I using 
the conservation equations 



It should be noted that there is no need for any 
arbitrary parameters 1 and 3. This is because the above 
dynamical equations can be determined exactly. 
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With Eqs. (9-58) and (9-59), the dynamic Eqs. (9-5^) 
and (9-55) and the grov/th Eqs. (9-51) and (9-52) become, 

2 



Kp + 



ll ^ R R ” ~T KP<T^> 



8 R "T Kp' 

2 



3 k_ ^ _ ML 

^ R^ R2 r2 



( 9 - 60 ) 



1 



<0,T) > + I - F 



1 k 



<T^> = - 



R2 



2R" R' 



R R2 



( 9 - 61 ) 






R" ^ 2 1 



R' 



R 



9p<T^> + <T^> 

K Rt 



2 2 
^-<0„T) > 



R 



( 9 - 62 ) 



3R<(aoT ) > + ^ <(90"^^ > " 



K ^ 2 

- 2 p9r<t^> 



(9-63) 



3 . Solution of the Moment Equations vjhlch Govern the 
Growth of the Fluctuations 

The growth equations are solved in Appendix K. 

The interested reader might note the differences betv;een 
the growth equations for the pressure free model and the 
radiative pressure model. The most significant difference 
is that the curvature constant k is in general not negligi- 
ble in the region where the matter density exceeds the 
radiation density. Therefore, in order to avoid difficult 
integration problems, the growth equations are solved only 
for the k equal zero universe. The advantage for doing 
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this is that the results can be expressed in terms of simple 
functions. (It seems that such problems presented by a 
nonvanishing k are inherent to nonuniform models. Sachs 
and V/olfe [2^1], for example, are also led to consider only 
Euclidean models for the same reason.) Nevertheless, the 
results for k equal to zero should have some relevance to 
k equal to plus and minus one universes. This is because 
the dynamics of the universe are nov; only slightly modified 
by curvature effects [25] and the solutions to the growth 
equations ought to therefore be approximately correct for 
such models also. 

With k equal to zero, the fluctuation free solutions 

2 

for Kp and R’ given in Appendix E are 



and 

R'" = ^ (9-6H) 

n 



where a^ is an integration constant different from a^ v;hlch 
appears in the analogous equations for the radiation filled 
universe . 

Then the growth equations become 



''' 2 “ 3<t^> 



'2 ' 



R^<OqT) > 



(9-65) 



and 



2 h 2 

3j^<OoT) > + r " 



f3a. 



— 3r<t^> 
r3 R 



( 9 - 66 ) 
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The solutions to these equations are given in Appen- 



dix K. 



- 7/2 -6 
— — = IqR " + + DqR 6 



<T^> = - I F^R-1/2 + I DqR-3 



(9-67) 

( 9 - 68 ) 



^ . Special Submodels 

Again, as in the radiation model, only the follov;ing 
three models are considered. 



(1) <T^> = i^R 



<Onf) > 



a 



= IqR 



-1 



(2 ) <T^ > 

(3) <T^> 



2 -1/2 ^ -7/2 

- ^ P^R : = F^R 



3 0' 



a. 



0 



-3 > _6 

^ D„R ^ ; = D^R ° 



9 0' 



0 



(9-69) 

(9-70) 

(9-71) 



The same difficulty occurs with the PQ-model as 
v;ith the the previous section. That is, it 

is impossible to find a nonzero Fg with a sign that is 
compatible with both of the Eqs. (9-70) (a^ is alv;ays greater 
than zero). Hence, only the Ig and Dg models are considered. 
It is to be noted that the Ig-model is a model for which 
the nonuniformity Increases with increasing R and the Dg- 
model is a model for v;hich it decreases with Increasing R. 

If Eqs. (9-69) and (9-71) are substituted into 
Eqs. (9-60) and (9-61) the follov;ing equations are obtained. 
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(Iq - fluctuations) 



Kp + 57 B(Iq)R 



-1 



3R 



R' 



3 B(Iq)R 



-1 



2R” R' 



R R^ 



(9-72) 



(Dq - fluctuations) 



tcp + 22 B(Dq)R 



-6 



3R* 

R^ 



where 



and 



- 2 B(Dq)R 



-6 



2R" R' 



R R" 









D„a 

B(Dq) = 



(9-73) 



99 - 7 ^) 



be 



The solutions of Eqs. (9-72) and (9-73) turn out to 



(Iq - fluctuations) 

3 a , 



Kp = - 



R' 



105 ^^^0^ 
2 R 



(9-75) 






(9-76) 



(Dq - fluctuations) 

3a , 



Kp = - 



20 B(Dq) 



R‘ 



a 



R< = ^ _ 



d 

R 



R 

2 B(D_) 



0 L 

R^ 



(9-77) 

(9-78) 
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If the same classification scheme, v/hich is depicted 
in Figs. 1 and 2, is used to classify these solutions that 
is used to classify the solutions for the radiation model, 
then it is easily found that the Iq - fluctuations are of 
type K and the - fluctuations are of type M. The discus- 
sion of this model and the radiation model is given in the 
next chapter. 
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X. SUMMARY AND CONCLUSIONS 



A. VALIDATION OP PROCEDURE 

Before discussing the results of this study regarding 
the development of a cosmological model to demonstrate the 
applicability of a moment hierarchy scheme to perturbation 
problems of general relativity, it vjould be well to show 
that the solutions obtained are compatible vjith ones obtained 
by a completely different procedure. 

Well knovm in the literature is the work by Sachs and 
Wolfe [24] concerning the anisotropy of the mlcrov;ave 3°K 
background radiation due to perturbations of geometry. In 
this study a scheme is presented for analyzing small, linear 
perturbations away from an unperturbed, k equal zero 
Robert son-Walker metric. The method of solving for the 
perturbations is to first take the spatial Fourier transform 
of the perturbations, solve them in k space, and then trans- 
form them back into position space. They do this for both 
the dust and radiation models. 

For the dust model (pressure free) they find two kinds 
of terms. In both types, 6p decreases although 6p/p increases 
for one and decreases for the other (in their notation, p 
is the mass density and 6p is its small perturbation) . 
Furthermore, the characteristic time for their relative 
increasing solution is of the same order as the character- 
istic time for the background metric. They report their 
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solutions for 6p to be of the form 



6p = DR + IR 2 (10-1) 

where D and I are functions only of the spatial coordinates. 
Because p goes as I/R , 

= DR"3/2 + JR (10-2) 

where the constant of proportionality Is absorbed Into D 
and I . 

As shown by Eq. (9-57), the ratio of 6p/p must be propor- 
tional to T. Hence, according to Eq. (10-2), t must also 
have the form 

T = DR"^/2 ^ (10-3) 

Thus, must be given by 

= D^R“^ + (DI)R“^/2 + j2r 2 ^ (10-iJ) 

But Eq. (9-68) says that 

<t2> = I Dq R"3 + R"^/2 ^ j^r2 _ (10-5) 

It Is seen that the exponents of each term agree vjlth 
those of Sachs and V/olf e . This agreement Is noteworthy 
because Eqs . (10-M and (10-5) are produced by two very 
different procedures. 

The determination of the characteristic time Tj for the 

Iq (Increasing) perturbation term t can be obtained from 

2 

the ratio of <t^> to <( 9 ^ 1 ) >. That Is, 



8l 



2 



1 



( 10 - 6 ) 



< T^> 

2 

<(a jT ) > 



„here c Is the speed of light. Using Eq. (9-69) 



R- 



(10-7) 






t- time T for the background geometry can 

The character istxc tim ig 

.. p p 2 (R') • That iSj 
be found from the ratio of R to ; 






1 _ rL- 

R'^ 



(10-8) 



With Eq. (9-69), Tb' 



7 . _ 



B 



R^ 

c^a. 



(10-9) 



indeed, the natic 6p/P which is pnopontional to t is pno- 
dicted py the »odel to ta.e place on the sa^e time sea 
as the background because the ratio of to Tg is o 

.his constitutes a second point oP agnee.ent with Sachs 

and Wolfe. models. 

are made only for tne uuo 

The above comparisons 

• „ between radiation models can be made an a 

But a comparison between that they neglect 

„ aachs and Wolfe indicate that tney 
3 imilar manner. c 

transport processes 

ponding model of this work is the 

which X and 6 are both equal to one. 

<^achs and Wolfe report that 
For their radiation model, Sachs 

(10-10) 

-fi 

6o = DR 
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so that 




2 -i4 
DR ^ 



( 10 - 11 ) 



where D Is a function only of the spatial coordinates. 

The general solution for <t^> given by Eq. (9-36) with 
X and 3 equal to one is 



<T 



V ^ R-V 



(10-12) 



The P 3_/2 solution is omitted from the general solution of 
Eq. (10-12) because, as stated in the text, certain sign 
problems are associated with it. 

But Eq . (9-2^) indicates that 



<Td> 



= - 2<t^> 



(10-lM 



This indicates that 
2 



That is. 



'd' 

.P. 



d 

IPJ 



a T ■ 



a 



'1/3 



+ 



D 



1/3 



R 






(10-15) 



(10-16) 



The comparison betvjeen Eqs. (10-11) and (10-16) indicates 
that although there is disagreement concerning the existence 
of an R term, there is agreement concerning the R term. 

It is not possible to assess the reason for this discrepancy 
owing to the omission of the solution procedure from Sach's 
and Wolfe’s paper. However, it is significant to this v.'ork 
that v/ith all the possibilities for the values of the 
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exponents due to the many possibilities for \ and g, there 
is an agreement concerning the existence of an R""^* term. 

Thus, there is a major agreement betv/een the functional 
forms of both the k equal to zero, dust and radiation filled 
models. Because this agreement is accomplished by tv/o very 
different perturbation techniques, support is given to the 
method of this work for treating perturbations by a hierar- 
chy of moment equations. 

It should again be mentioned that the agreement of the 
radiation model is contingent upon the selection of the 
"zero leakage" values for X and B because Sachs and V/olfe 
assume that the velocities are zero. It is reasonable to 
ask how valid is this assumption. In the derivation of the 
model, the attempt was made to include the effects of the 
velocities by the- introduction of the parameters x and g 
whose magnitudes are estimated in Appendix H to be of 
order one . 

Figures 3 and ^ give the types of solutions that can be 
expected for all possible combinations of X and 3 . In the 
case of the I 2./3 “ perturbation, the model classifica- 

tion for X and 3 values of order one is clearly type K. But 
for the “ ^^P®’ both the K and M - type solutions are 

possible. In fact, the zero leakage values of X and 3 are 
observed to lie on the boundary betv/een the K and M - type 
solutions. Thus, if a " '^^P® >iiodel is being considered 

the neglect of the velocities is not trivial for very differ 
ent dynamical situations could occur even if the velocities 
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are extremely small. But for an l^/S ~ neglect of 

the velocities does not change the solution type. In this 
case the approximation of zero velocities seems to be justi- 
fied. Indeed, Sachs and V/olfe consider only this type. 

However, A and B not only affect the classification but 
also determine the exponent (v;hich is j + ^< for and i+4 

for the expression for <t^>. If Figure 5 

the lines of constant i and j are given in A-B space. Prom 
these curves, it appears that for A and B of order one, 

-3 < j + $ - 8 

3 :S i + 4 < 8 . (10-17) 

In the case of A and B equal to exactly one, j+^ is equal to 
-4 and i+i| equal to These estimates are conservative. 

Perhaps a more reasonable estimate would be to assume that 
the likelihood of A and B both being on the large side of 
their order one magnitude is small. Then the upper limit 
might be pushed down to perhaps -6 and 6 respectively. 

Thus it seems that the zero velocity assumption can 
reasonably be questioned on the grounds that the exponent of 
<T^> could be different enough from the zero velocity expo- 
nent that a many order of magnitude change of R could produce 
a value of <t^> that is much larger (or smaller) than the 
corresponding zero velocity value. This can be seen more 
quantitatively as follows. Consider, for example, a radia- 
tion model with an perturbation. At some value 

of R equal to R, 
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(10-18) 



<T^> 



= 1/3 



r1 + ^ 



where the sub and superscript i refers to the i-value of the 
exponent. If it is assumed that the ratio <t^>^/<t^>q is 
one at R, then this ratio is in general different at other 
values of R. That is, this ratio is a measure of the effect 
of a given combination of A and 3 upon the nonuniformity 
in terms of the "zero velocity" nonuniformity. This ratio 
is easily found to be, for an arbitrary R 



<T^>. 



<T^>, 



R 



1 



R 



(10-19) 



where i = - ^ + 2 /2 ( A + 3 ) . 

Thus, unless A and 3 are such that i is near enough to 
zero so that <t^> is of order the effects of the 

expansion through many orders of magnitude upon the nonuni- 
formity compared to the zero velocity nonuniformity is signif- 
icantly increased as seen from Eq. (10-19). Ifj hov;ever, 

R changes at most by only a few orders of magnitude in a 
given problem, the discrepency is not as serious and may be 
tolerable if i is close enough to 0. Hence, although the 
zero velocity approximation is extremely convenient mathe- 
matically, its validity does not appear to be on solid footing. 



B. DESCRIPTION AND COMPARISON OF MODELS 

As stated in Section A of this chapter, there are two 
intrinsically different types of nonuniformities associated 
with both the dust and radiation filled models. One is an 
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increasing (I-type) and the other is a decreasing (D-type) 

2 

nonuniformity. If the ratio of <x^> to <(8 jT) > is computed 
for all four models (tv;o for each a value) in the same v;ay 
as done previously so that a characteristic time might be 
estimated for the relative perturbations, then it is found 
that, for all models, this characteristic time is of the 
same order as the characteristic time for the average 
background . 

Because of the existence of the two D and I-type pertur- 
bations, it is natural to ask which kind characterizes our 
universe. This question can only be answered through obser- 
vation. Although such observations are difficult to perform, 
there seems to be some suggestion that the I-type of pertur- 
bation is the one particular to our own universe. One 
reason why this might be true is that the average background 
metric might be expanding more rapidly than the supercluster 
to which our galaxy belongs. A second reason is that, in the 
case of D-type perturbations, galaxies might be discarded 
from their otherwise clustered configurations [2^]. 

Granting the I-type dust model, it is possible to estimate 
the parameter Iq which governs the amplitude of the metric 
perturbation <t^>. This can be accomplished by using the 
data reported by Partridge and VJilkinson [19] on the aniso- 
tropy of the so-called cosmic blackbody radiation. They 

-7 

report that the ratio 6T/T is of the order 10 where T is 
the temperature of the radiation (3**K) and 6T is the fluctua- 
tion of T with angle of observation. Sachs and V/olfe [2^] 
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give a very simple expression that relates 6 T/T to 6p^/p^ 
where Pj, is the current mass density. 



6 T 

T 







( 10 - 20 ) 



Ro’ 

where H_ - 5 — (* denotes derivative v.Tt . x°), R is the 

n n Q 0 

current value of R, and L is the characteristic length of the 
perturbations. By Eq . (9-57), Eq. (10-20) can be written as 



_2 _ 1^ 
9 



(HpL)' 



6 T 



The constant can be estimated by means of Eqs 



and (10-21) . 


That 


is. 




I^ - 


e" 


'6t]' 


1 

fl6] 


0 


Ro" 


tJ 


1 9j 



(H^Ro)^ 



(10-21) 

( 9 - 69 ) 



( 10 - 22 ) 



where 




(10-23) 



If H„R is taken to be one and if the factor 16/9 is neglect- 
n 0 

ed since it is of order one, Eq. (10-22) becomes simply, 



I 



0 ~ 




( 10 - 2 ^ 1 ) 



If the fluctuations are assumed to be the size of about 
1000 MFC [2^, 31 ] (one MFC is equal to 3*26 light years) 
so that c is of order ten (R^^ ~ 10^ cm), then for Fartridge 
and Wilkinson's estimate of 6T/T, Eq. (10-2^1) gives for Iq 
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(10-25) 



I 



0 



10 



-58 



-2 

cm 



Using this value for Iq, it is interesting to determine 
the effect a term of this strength v/ould have upon the average 
dynamics if it is assumed that the long v/avelengths are 
important to the expansion characteristics [35]. It is v/ell 
known that the k = 0, unperturbed uniform universe model 
expands forever [2]. (This can be seen easily by integrating 
Eq. (9-6^) to find R a But ov;ing to the additional 

source term, the dynamics of the model is perturbed. As 
pointed out earlier, the solution for the k = 0, Iq model 
is type K. Because this means that R is perturbed below 
the uniform solution, (see Fig. 1), the question naturally 
arises v;hether or not there exists a maximum value for R. 

To answer this, R' must be set to zero in Eq. (9-76) and a 
positive value of R greater than R(, must be determined. If 
this is done, such a value R is found such that 

^ m Y 



max 



16 

31, 



1/2 



Using Eq. (10-25) , 



max 



10 cm . 



( 10 - 26 ) 



(10-27) 



That is, if the observed anisotropy of the blackbody radia- 
tion is due to perturbations of space as large as super- 
cluster distances and are of the Iq - type, then the 
"effective energy" of these longv.^ave fluctuations is 
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sufficient to prevent an unlimited expansion of a k = 0 , 
dust filled universe which v;ould otherwise expand forever. 

It should be pointed out that this model is not valid 
for all values of R. For small R, the radiation density, 

which goes as R , overv/helms the mass density, v;hich goes 

-3 

as R . At this point, it is no longer valid to speak of a 
dust model and it is proper to use the radiation model for 
smaller R-values [15]. On the other hand, for large R, 

<T^> approaches one. As it does so, the assumption of small 
amplitude fluctuations which was used to truncate the moment 
hierarchy, becomes invalid. The point at v;hich this happens 
can be easily determined. V/ith <x^> equal to one, Eq. (9-69) 
gives 



R = 
c 






1/2 



(10-28) 



where R is that value of R for v;hlch <t^> is equal to one. 
c 

If Eq . (10-28) is compared to Eq . (10-26), it is found 

that R <R which indicates that at the "turning point," 
c max 

the fluctuations are as lai’ge or larger than the average 
metric. However, it should be noted that because the model 
is not applicable at the turning point itself, it may not 
therefore be justifiably used to determine This diffi- 

culty is not serious though since the ratio of the "extrapo- 
lated value" of R given by Eq. (10-26) to R is 

ITla.X ^ 



^max ^ i (10-29) 

^c /3 
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should 



which is of order one. Thus the real value of R, 
still be of the order of 10^® cm as given by Eq. (10-27) 
since the theory breaks down close to the estimated value 
of R .It should finally be noted that if other models 
are considered v/hlch also have "turning points," (either as 
a minimum or maximum) the v/eak fluctuation assumption breaks 
down just before the turning point occurs. Thus, if the 
physics of an assumed turning is to be accurately discussed, 
a strong fluctuation model should be used. 
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APPENDIX A 



The Reynold Stresses 

The Navier-Stokes equation for an incompressible fluid 
is [12] 

ii + v.vv = -^vp + vV^V (A-1) 

where V and P are the velocity and pressure fields relative 
to a fixed coordinate system and v is the ratio n/p where 
n is the viscosity coefficient and p is the mass density. 

For Incompressible flov;, p is a constant throughout all 
space . 

In order to gain some insight into the perturbation 
technique that is presented in this paper, the reader should 
look upon the Navier-Stokes equation as if it v;ere the 
Einstein Field equations. The left hand side of Eq. (A-1) 
is analogous to the geometrical side of the field equations 
and the right hand side is analogous to the mass side 
(that is, the side which contains the energy momentum ten- 
sor). The velocity field V plays the same role in the 
Navier-Stokes equation as does the metric tensor g^^ in the 
Einstein field equations. 

In fluid turbulence theory, the following decompositions 
are assumed [^]. (In this appendix-, the equation after the 
semicolan refers to the analogous equation in the text.) 

V = <V> + V (A-2) ; (^-9) 
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f 



I 




p = <p> + p 



(A-3);(iJ-9) 



where < > denotes "average." 

If Eqs. (A-2) and (A-3) are substituted into Eq. (A-1), 
the following equation results 



+ <V> • V<V> + V • V<V> 



+ <V> • Vv + V • Vv 



11 ^ 

= - — V<P> - — Vp + vV^<V> + vV^v 
P P 



(A-M 



In fluid turbulence theory it is assumed [^] that if 
A is a random field, 



and 



' It 



<iJ^ = <A> 



(A-5);(B-2) 



(A-6);(B-2) 



3x‘ 



9x' 



Because of Eqs. (A-2) and (A-3), 



<v> = 0 



and 



<p> = 0 



(A-7);(^-l^) 

(A-8) 



Then, Eq. (A-^) becomes 

+ <v> • v<v> 



3<V> 

at 



- - V<P> 
p 



+ vV^<V> - <v • Vv> . 



(A-9);(6-8),(6-9) 



It should be noted that Eq. (A-9) is symbolically iden- 
tical with Eq. (A-1) with V and P replaced with <V> and <P>. 
The only difference, hovrever, is the appearance of the extra 
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source term from the "velocity side" v;hich results from the 
nonlinear V • VV term. It is, in general, nonvanishing. 

This extra term is called the Reynold Stress. It can 
be interpreted as an additional source term (in addition 
to the pressure and viscous sources) which acts on the 
average flow field (<V>). Such an Interpretation can be 
made because the left hand side of Eq. (A-9) is symbolically 
identical with the left hand side of Eq. (A-1). 

In order to obtain the functional form of the Reynold 
Stresses, a hierarchy of equations can be built up from the 
Navler-Stokes equation of Eq. (A-1) by multiplying it v;ith 
various orders of V and averaging. But the nonlinear term, 
as already seen, alv;ays produces a moment of one order 
greater than the order of the moments which comes from the 
other linear terms. Thus, there is always one more unknovm 
than there are necessary equations available to solve the 
system of equations. This problem is the well known closure 
problem of turbulence theory which is discussed in the 
text [3]- 
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APPENDIX B 



Definition of the Averaging Operation 



There are various definitions of average v.'hich are pos- 
sible. Many theorists v:ho v/ork in the field of fluid turbu- 
lence use an integral expression for their definition of 
average. For example, Curie and Davies [4] use a time 
integral. This is an acceptable definition for experimental 
work but as pointed out by Lumley and Panofsky [131 j such 
an integral expression is not always easy to define in a 
meaningful v;ay. They suggest that an ensemble average is 
probably the most basic. Thus, the notion of an ensemble 
average is taken as the definition for the averaging process 
< > for this model. 

Suppose D(x^,^) is a stochastic process v/here x^ are the 
four space-time coordinates and ^ is the ensemble parameter 
which designates D as a possible outcome for an experiment 
(see also Ch. IV). Then <D> is defined as 



<D> 



lim 



D(x^; Cj) + ... + D(x^; 
_ 



(B-1) 



where <D> is only a function of x'^ and is not dependent on 
5 . This is the same definition as that used by Lumley and 
Panofsky . 

Because this definition involves a summation, the follow- 
ing equation is assumed valid. 
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3 <D> 

y 



<3 D> 

y 



(B-2) 



Suppose it Is desired to know the probability of finding 
D(x^) falling below a certain value <S . This can be deter- 
mined by counting the proportion of members of the ensemble 
for which D(x^) 5 6. From this a function P can be defined 
as follows 

F(6;x^) = P{D(x^) 1 6} (B-3) 

where P stands for "Probability that." The probability 
density function is defined by 

„ 3F(6;x^) 

f(D;x^) = . (B-^) 

Then the average of any function h[D(x^)] can be written as 

00 

<h[D(x^)]> = / h(6)fU,x^)d6 . (B-5) 
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APPENDIX C 



An Implication of Spatial Stationarity of the Random Field 



An implication of statistical spatial stationarity is 
rather easy to state because of the assumed form of the line 
element. On a hypersurface of constant x° , dx° is zero so 
that by Eqs . ( 4 - 5 ) and ( 4 - 6 ) 



ds . = - 






do‘ 



1 + 



k r 



2 1 



(c-1) 



It is pointed out in Chapter IV that Rj^ has the meaning of 
"local radius of curvature." It is also a stochastic process 
which means that at each point x^ of space time it has a 
probability density function associated v/ith it (see Appen- 
dix B) . 

In order to bring about a simplification of the proposed 
model, consider the follov/ing assumption. Although the 
hyperspaces of constant x® are clearly nonuniform, and under 
go various expansions at each point, it is assumed that the 
probability of finding a given expansion at any point is 
independent of position on that hyperspace. This require- 
ment is then one of statistical spatial stationarity, a 
concept v;hich is familiar in fluid turbulence theory. This 
implies that the probability of measuring a given value of 
Rj^^ between say R^^^^ and Rj^^ + dR^^^^ ought to be independent 
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of the coordinates on this spatial hypersurface of constant 
time x° because this random function contains all the expan- 
sion characteristics. Equivalently, the probability of 
finding t between T and T + dT ought to also be Independent 
of the spatial coordinates. Hence, the probability density 
function f(T;x^) defined by Eq. (B-^) must be such that, 

f (T;x® ,x^^ ) = f (T;x° ,x*^ ' ) 
i i ^ 

for every x'^andx where j = 1, 2, 3- This can only mean 
that [l8] 

f(T;x^) = f(T;x°) (C-2) 

That is, f can only be a function of the coordinate x® and 
T. But Eq. (B-5) says that for any function h of the process 

00 

<h(x)> = / h(T) f(T;x°)dT . (C-3) 

— 00 

Hence, <h(t)> can only depend upon x° and not upon r, 6, (P . 

It follows that 

9.<h(x)>=0 (C-^) 

J 

where j = r, 6, . 
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APPENDIX D 



The Form of the Energy-Momentum Tensor 

The general form of the energy-momentum tensor for a per 
feet fluid is given by [1, 29] 



density and pressure. That is, these are the density and 
pressure vjhich vjould be measured by an observer moving 
instantaneously v/ith the fluid. Tolman [29] defines a per- 
fect fluid to be a mechanical medium incapable of exerting 
transverse stresses in the proper frame. From this defini- 
tion follows Eq. (D-1), 

Tolman [28] shows that the energy-momentum tensor for a 
disordered electromagnetic field, v;hich is expressed in term 
of the electric and magnetic field strengths, can be viewed 
as a perfect fluid for v/hich 



However, if the pressure is due to the random motion of 
particles of cold dust, then the equation of state, which 
is well-known to kinetic theory, is given by [1] 



= DUkU^ + ~ - g^^) 



(D-1) 



c 



where = dx^/ds and D and P are the proper hydrostatic 



P = I c^D. 



(D-2) 




(D-3) 
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where is the root-mean-square random velocity of the dust 
particles. But this velocity is much less than c so the 
pressure, for all practical purposes, is essentially zero. 

The two models, therefore, can be written together in 
one equation. 

P = ac^D (D-4) 

where a = 1/3 if universe is radiation dominated 
a = 0 if universe is matter dominated . 



Now, the expression for the mixed energy momentum tensor 
is needed. This is easily found by means of the metric 
tensor g 

V 



= DU^U + — (U^U 

V V ^ 2 V 



where 



U 

V 



g U 
^av 



a 






(D-5) 



(D-6) 



In order to evaluate Eq. (D-5), and need to be 
determined. This can be done as follov/s. Consider a fluid 
element which moves along the trajectory x'^(x°). Along 
such a trajectory. 



ds^ = (dx“) 



1 - e 



M da 
(dx®) 



With the definition 



u = e 



M/2 



dx‘ 



(D-7) 



(D-8) 
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equation (D-7) becomes 



= + (1 - . (D-9) 

dx° 



Notice that if the motion of the fluid element is along a 
null geodesic, u^ must be equal to one. In order to keep 
the interval of such trajectories timelike (interval real) 
so that the local velocity of the fluid element does not 
exceed the light velocity, it is obvious that from Eq. (D-9), 

0 < u2 < 1 . (D-10) 



It is natural to identify (uc) with the local velocity of 

the fluid element at the point x^ because e da is the 

physical distance that the fluid element moves in a time 

— dx ° . Now define 
c 



Y 




(D-11) 



SO that 

uo = = Y . (D-12) 

ds ' 

Therefore , 

uJ = (D-13) 

ds dx° ds dx° 



But Eq. (D-8) can be written as 



u^ = e 



M 

e 




fdr 1 


2 

+ r^ 


fde 1 


+ r^sin^e 


rd(j, 1 


2 ' 






dx° 




dx 




dx® 

k J 


- 



. (D-lil) 



If the quantities u^, u^ and u^ are defined such that 
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_ M/2 dr 

® “ 5 " 

^ dx 



u = re 



M/2 d6 



dx 



then 



. . M/2 d(j) 

u = r sin0 e — i— 

4> dx“ 



u" = u ' + U.2 + u 2 



6 (p 



It is natural to Identify as the j 
local velocity at x*^ . Hence, 



th 



(D-15) 

(D-16) 

component of the 



U® = Y 

„i ^ -M/2 
= Ye u 



-M/2 



U* = Y - 



-M/2 



u 



r sin6 4> 



(D-17) 



By Eq. (D-6) , 



U = Y 
0 



U = 
1 



V M/2 
- Ye u 



U = 
2 



V M/2 

- Yre u. 



U = 

3 



M/2 

- yr sin6 e ^ u 



(D-18) 
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Then Eq. (D-5) becomes. 



T ^ 

0 


= D + — u2 
c ^ 


Ti 

1 


= - Du ^ (1 + u " - u2) 

^ c2 ^ 



= - Du^- ^ (1+ u/- u2) 

0 C 2 6 



T ^ 

3 


= -Du/- — (l + u^-u2) 

4> C 2 <P 


1 


= - = - (D . 

C2 ^ 


T® 

2 


= - r2e^T^ = _ (D + ^) re^^^u. 

0 ^2 6 




= - r^sin^e e^"^T^ = - (D + — ) r sin6 e^^^u, 

» c2 




_ = _ (D + Z_) u. 

1 ^ 0 


1 

T 

3 


= r^sin^S T^, = - (D + — ) r sln0 u u, 

1 ^ ‘f’ 




= sin^0 T^, = - (D + — ) sin0 u„u, (D-19) 

c2 0 <*> 



where D = 

P = PY^ . 

To produce Eqs.(D-19), the follov.'ing identities are needed. 

- 1 = 

Y^u! + 1 = Y^(l + u - u^) . (D-20) 

u J 
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APPENDIX E 



Solution of Cosmological Equations v.^ith a l/R*^ Source Term 



The tv;o usual cosmological equations of general rela- 
tivity with zero cosmological constant but v/ith a 1/R^ 
source term are [1] 



and 



where 



Kp + Kp^ = 



KP 



3k _ 3 ^ 

R2 



R2 

k 2^' _ RJ_ 

R2 R r2 



dx' 



(E-1) 



(E-2) 



p 

V ^ra 


(E-3) 


P = ac2p ■ 


(E-4) 


P = a c ^p 

V V V 


(E-5) 


8itG 

K = - 


(E-6) 


c^ 




. = 


(E-7) 



The difference between Eqs . (E-1) and (E-2) is 









R 



(E-8) 



The derivative of Eq. (E-1) is 

6kR' 



K(p' . p;) = 



- 6 



fR’l 


fR'l 


R 

V 


[r J 



(E-9) 



1QI\ 



The substitution of Eq. (E-8) into Eq . (E-9) gives 



(p’+p„ ) + (p+p„) + ^ — = 0 



V' R 2 

C ^ 



If Eq. (E-8) is multiplied by R^, it becomes 



((P+P^)R') + 



(P+P„) 



V 



(RM = 0 . 



(E-10) 



(E-11) 



This is the conservation of mass equation v^here and 

are the mass density and pressure of the source respectively. 

It is valid for any equation of state P(p) and P^(p^). 

But with Eq. (E-5), P^ = a^c^p^ and Eq. (E-4), P = ac^p 
Eq. (E-11) becomes 

(pRM'. + ap(RM* = - (p^RM' - a^p^(R'M' . (E-12) 



The homogeneous equation associated v;ith this equation is 

(pR^) + ap(R^ = 0, (E-13) 

the solution of which is 



P 



A 

j^3(a+l) 



(E-1'4) 



In order to solve Eq. (E-IP), Eq. (E-1^) may be substi- 
tuted into it assuming A is now a function of x® . Eq. (E-12) 
becomes 



A' = - p^’ - 3R^“'^^R'P^(1-Hx^) 



With p = — , Eq. (E-15) becomes 

V j^m 



(E-15) 
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1 ! 




fm-3(l+a )'l N ' 

V r R 3 (a + 1 ) -m'l 

T^-3 1+rvU ) • 



(m-3(l+a)j 

But Eq. (E-16) Is easily integrable. It turns out that 



(E-16) 



where 



A = 



C = 



- + B 



m-3(l+a^) 



m-3(l+a) 



(E-17) 



(E-18) 



Finally, the substitution of Eq. (E-17) into (E-1^) gives 



B _ CV_ 

^ " j^3(a+l) pjOi 



(E-19) 



The substitution of Eq. (E-19) into Eq . (E-1) gives an 



expression for R' . 



R' 



K 

3 



B V(l-C) 

p3o+l j,m-2 



- k 



Thus , 



and 



Kp = - 



3a 



p3(a+l) p^rn 



R 



a ^ b 



where 



a = - 



R3a+1 pm-2 



KB 

3 



- k 



f = KCV 



^ _ kV(C-I) 
^ - “1 



(E-20) 

(E-21) 

(E-22) 

(E-23) 

(E-2i|) 

(E-25) 



If Y = 0 (no source term), then the solutions become 
the usual cosmological solutions for p and R. For a radiation 
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filled universe (a = 1/3) [29] 



3 



and 



Kp = 



R- 




(E-26) 



(E-27) 



For a pressure free universe (a = 0), the zero source 
term solutions are 



and 








(E-28) 



(E-29) 



It should be noted that the meaning of the constant a 
is different In the case of the radiation model (denoted 
by a^) than In the case of the pressure free model (denoted 
by ) . 
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APPENDIX F 



Elimination of the Spatial Derivatives 
from the Dynamic Equations for a Radiation Pilled Universe 



The terms involving the spatial partial derivatives of 
T can be removed from Eqs . (9-1) and (9-2) by the substitu- 
tion of Eq. (9-3). The result is 

“ 2 ^ + ~ <T^> 

R ^ 



= + 1^’ 12Rji 

R2 R r2 

- “ 2 i~ I" ^ ^ ^ <T^> 



iiii + 3R_L' ^ ML 

R^ R R^ 



(F-2) 



But the right hand side of Eqs. (F-1) and (F-2) are not 
the expressions that would result from substitution of the 
average metric into the Einstein field tensor. In order to 
bring these equations into such a form, Eq . (F-2) must first 
be multiplied by minus five and added to Eq. (F-1). If this 
sum is multiplied by three eighths, then the desired form 
results . 






3 k 3k 3R* 

“ ^ R^ 
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(F-3) 



To obtain the other equation, Eq . (F-1) must be multiplied 
by minus 5/3 and Eq. (F-2) by 3- IF these tvjo equations 
are added and divided by eight, the resulting equation is 
formed . 



K£ 

3 



3 

TE 






1 



^ a 
R ^ 0 



<T^> 



1 



<( 3 



2 

> 



1 k 

¥ ~ 
R- 



<T^> 



k ^ 

R2 R r2 



2 

The term 3 <x^> can be replaced by using Eq. (9-^)- 

final pair of equations are 



(F-iO 

The 



5 

Kp + ^ K<xd> + 






> + 



3 

F 



R 



3„<t^> 



3 k. 
^ r2 



<T^> 



3k _ ^ 
R2 r2 



- f - 1 «Td> 




<(3^t) > 



1 R 

F 



3o<T^> 



1 k 
^ R^ 



<T^> 



k 2^' _ R/ 

R2 R r2 



(F-5) 



(F-6) 
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APPENDIX G 



Elimination of the Spatial Derivatives 
from the Dynamic Equations foi"" a Dust Pilled Universe 



The elimination of the terms which contain the spatial 
partial derivatives in Eqs. (9-^8), (9-^9) and (9-50) 
gives 

- f Kp - I I” ^ " 2 

+ ^<,2, = 12k ^ 1^ ^ 1^- 

r2 R 



5 1 5_* ^ 

6 P 2 R 



<t'> - I 3o<t^> + ^ <Oo't)^> 



I- K<Td> + — 
^ R2 



<T^> = 



iiii + iiLL 

R^ 



R' 



3R*' 

R 



(G-2) 



Now, if Eq . (G-2) is multiplied by minus five and added to 
Eq. (G-1) and the resulting equation multiplied by three 
eighths, the following equation is obtained. 



Kp 






+ 



5 

F 



K<xd> 



3 k 



<T^> 



3k _ 3^ 

R^ R^ 



(G-3) 



If Eq. (G-1) is multiplied by minus five thirds and 
Eq. (G-2) by three, and if these tv;o resulting equations 
are added and divided by eight, the following equation is 
obtained . 
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- li + 



1 R' 

? r -0 



8 <T^> 



1 2 c: 

^ <(8 .t) > + ^ 






K<xd> 



1 k 

W ~ 

R‘ 



<-f 2 > _ 



R^ 



2R" R' 



(G-lJ) 



R R = 



Now, the term 8^<T^>that appears in Eqs. (G-3) and (G-4) 
can be eliminated by using Eq. (9-5D* The result is 

Kp + ^ <(8 ijT) > + ^ ^ ¥ <<Td> 



R^ 



3k _ 3RJL 

r2 R^ 



(G-5) 






Ik , 2, _ 

IT — > = 

R2 



l£_ 

R2 



2R" _ RJ_ 
R R^ 



(G-6) 
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APPENDIX H 



The Dynamics of the Energy Momentum Tensor 
for Disordered Radiation 



The dynamics of the disordered radiation fluid is 
described by Eq. (9-10) 

T^^^; V = 0 

where 

T*^^; V = 9 T’^'^ + T^^“ + T°‘^ . (H-1) 

V va av 

The components of T are given in Appendix D, and the com- 
ponents of are given by Eq . (5-10). Because of the 

a 

summation which appears in Eq. (H-1) with respect to the 
indices v, Eq . (H-1) represents a set of four differential 

equations. Adler [1] points out that the equation v;ith y 
equal to zero is the energy-conservation equation v;hile 
the remaining three are the relativistic equivalent of Nev;- 
ton’s law of dynamics applied to an ideal fluid. The con- 
servation of energy is therefore expressed by 



g ipOV ^ pV ipo*^ + pO ijiOtV _ Q 

^ va av 



(H-2) 



With the appropriate substitutions from Appendix D and 
Eq . (5-10), Eq. (H-2) becomes 

8 (D + — u'n + i (9„M)[(D + — )(3 + u2)3 

0 9 Cl u J 



+ e-”/29^[(D + ^ ] 
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(H-3) 
cont . 



+ e ^^(9 M + |)[(D + — )u ] 



2 r ' 



cont . 
(H-3) 



c 



-M/2, 1 

r ''6- 



+ e i 9,M + £2 ^)[(d + L^, uj 



^ r-i^ 3 [(D + u ] 

c ^ 



+ e -”/2 ( 1 



r sinT “4,1 ' « 



If the equation of state given by Eq. (5-22) is used, 
Eq. (H-3) becomes. 



9^D + 29^MD -I- j 9 ^(DuM + | (9^M)Du2 



+ ^ [e ^'^^a^(Du^) -I- e ^^^^(9^M 4- -)(Du ) 

r 



+ e - 9n(Du„) -H e ( - 9,,M + ^ )(Du^) 



6 ' 0 



* rwT 5/%’ 



^ <FTnre 

In order to find suitable expressions that might replace 
the correlation functions <xd> and <x9j|d> in Eqs . (9-^) and 
(9-5), the following procedure is used. 

If Eq. (H-i|) is multiplied by a special function 
h(x,9gX) and averaged it becomes 

<h9oD> + 2<h(9^M)D> + j <h9g(Du^)> 
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(H-5) 



+ j <hDu^OoM)> + H<e 



r sin 8 ^(}) 



9a(Du, )> 

<P <P 




(9.M)> = 0 . 



It should be noted that in Eq. (H-5) only partial deriva- 
tives with respect to x° and c}) appear explicitly. Those 
terms which contain partial derivatives v?lth respect to r 
and 0 are buried in the terms containing the <p partial deri- 
vatives. To see how this can be done, consider the follov;- 
ing argument. In Eq. there are three types of terms 

involving spatial partial derivative v/ith respect to r, 0, 
and (}). On a spatial hypersurface of constant x® it is 
reasonable to assume that the average of a quantity associa- 
ted with a given direction in space ought to be equal to 
the average of that same quantity associated v;ith a different 
direction. Since the r, 0, and (J) terms of Eq . (H-^) have 
identical roles, (i.e. transport of mass in r, 0, (j> directions 
as Adler [1] points out) then the average of these terms 
(multiplied by the function h which has no explicit spatial 
dependence) should be equal. 

Now, D can be decomposed into its average and fluctuat- 
ing components. That is. 



such that <D> = p . 

If Eq. (H-5) is expressed in terms of x instead of M, 
it becomes, using Eq . (H-6) 



D = p(x“ ) + d(x° , r , 0 , (}) ; C) 



(H-6) 



Il4 



9oTT<h> + <h9pd> + p<h> + ^^<hd> 



R 



+ 2 p<h 






il + Tj 



> + 2<h 



U + Tj 



d> + i 9pp<hu^> 



+ ^ P<h9o(u^)> + I <hu^(9jd)> + ^ <hd(9jjU^)> 



4 - R 



2 -. 



+ -r- P 5 — <hU > + : 5 - P<hU 

on 0 






1 + T 



^ R' — 7 2 - • 

+ I <hdu^>* + j <hdu' 






,1 + T, 



+ :,<e-M/2 3 p> , Hp<e-"/2 __h , 

r sin 6 (j) r sin 6 (^ <{) 






Tv. , h -M/2 hd ^ 

n °Ad> + ^<e 3 T 

r sin 69 n sin 6 4 > 



+ ^P<e 



hu 



<P 



r sin 0 






1 + T 



-M/2 

+ /|<e ^ 



f \ 



V sin 6 



1 + T 



> = 0 . 



(H-7) 



It should be expected that weak fluctuations of geometry 
do not create pressure gradients which accelerate the radia- 
tion fluid to local velocities near the speed of light. 

That is, 

U^<<1 



whenever 

T^<<1 



(H-8) 



This makes certain moments of Eq. (H-7) small compared to 
others . 

Consider now tv;o cases. Let h symbolically stand for 
two special functions, x and 9^1. Then 
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<h> = 0 . 



(K-9) 



Thus, If the higher order moments of Eq . (H-7) are neglected 
by virtue of the smallness of u^,this equation becomes 

<h9jd> + — ^ ■■ <hd> + 2p<h9jT> + ^ p<h9oU^> 



+ I p<hu^9 T> + l4p<e ^ S.u,> 

3 0 r sin 6 <}) (J) 

hu 



+ 4<e 



-M/2 3 ^ -M/2 

r sin d (p 



4 > 



r sin 0 (}) 



9 ,T> 



= 0 . 



(H-10) 



As pointed out by Harrison and others [9], any perturba- 
tion of the radiation fluid characterized by d, u^, Ug, u^ 
will propagate with the sound velocity of the fluid. For 
radiation this velocity is near the speed of light. 



V 

sound 



= A 

/ dD ^ 



c 



Hence, such an accoustlcal signal v\’Ould propagate along 
"near" geodesic lines. That is, ds^ (accoustlcal signal) 
is near enough to zero that 

dx° e'^^da - r sin 6 d4> • (H-11) 



But this implies that the three moments v.’hich contain the 

spatial partial derivatives with respect to 4> in Eq. (H-10) 

-M/2 

ought to be of the same order as those for which e x 

^ Q 9c}; is replaced by 9(,. For example. 



-M/2 h 

r sin 0 



<p (f) 



<h3.u^> . 



Il6 



Then, Eq. (H-10) becomes after discarding the insignificant 



moments 



— ^R’ — — 

<ha d> + — 5 — <hd> + 2 p<ha T> 
on *^0 



. -M/2 

+ Mp<e 



r sin 0 (p (p 



9 u > = 0 . 



(H-12) 



If the first two terms of Eq. (H-12) are consolidated 
into one term fj<h9^d>(for an appropriate fj), it becomes 

f <h9 d> + V — ^7 

1 0 ^ ^ ^ ^ I , 



r sin 0 (})(}) 



= - 2p<h9^T> . 



(H- 13 ) 



If the characteristic time for the fluctuation is comparable 

1 R 

to the characteristic time — ^ of the average geometry, 
then the second term is equal in magnitude to the first. 

In this case, f j , which is a function of x° , is of order 
one . 

Nov; consider the remaining equations v = 0 for j 

equal to one, tv;o, and three. After some rather tedious 
calculations, these equations become 



V = 0) 



3 r(D + ^) + (| 9„M)[(D + e-’^^/^u ] 

0 o 1 ’ ^ V 2 



+ 3„{e-"[Du^ + ^(1 + - uU]} 



+ (I 3 M + 3) e-‘‘[Du^ + i_ (1 + - uU] 

2 r ^ c 



P 



- J (9 M + 2) e-^[Du^ + ^ (3 _ 2uM] (H-li^) 

^ r cont. 
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cont . 



-M 

+ p— 



.5 P e-” 

+ (| 3qM + cote)[(D + — ) § U^u^] 



r re 



D -M 

+ 3 [(D + — ) ^ ^ 

4) ^2 r sine r (p 



-M 



+ I (3,M)[(D + — ) - - -- ., u u ] = 0 
^ (j) ^2 ^ siriG r (p 



V = 0) 



3 J(D + — ) — u^] + (4 3 „M)[(D + — ) - 



r 6 2 0 



2 r 



-M 

+ 3 [(D + — ) ^ u u^] 

r ^2 r r0-^ 



c Zj P e ^ 

+ (?■ 3 M + ^)[(D + ^) ^ u^uj 



2 °r r- 
-M 



r r 0' 



+ ^ 3 



2 . P 



r e r 



[Du! + ^ (1 + u, - u2)] 



® c2 



t (| i 3e« - j ^ [Du^ + fr <1 " % 

“M 

- I (i 3gH) ] 2_ [Du* + (3 - 2u*)] 



cote 



-M 

^ [Du" + ^ (1 - u" - u2)] 



-M 

+ 3J(D + 2_ UgU^] 



r sine (j) 



(H-li^) 



- u^)] 



(H-15) 
cont . 
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If 

state 



5 1 



^2 r slnG + ~) V” ^0^6^ ~ ° (H-'is) 



-M 



cent . 



2 ^ r "e^(J)- 



(T^^; V = 0) 



P e~^/2 

aJCD + — ) ^ u^ 



. 5 P e-^/2 

] + (| S,M)[(D + ^) 

c 



p 

^ »rt(D 



f (| S^M + 1)C(D + ^) ir^u^u,] 



-M 



1 . r/^ . P X e 



-M 






— M 

+ (I ^ 8,M + 3- - q . gt9 )[(D + L.) u^u^] 



'2 r e 

__1 3 

r sin0 ({) 



+ (I ^ 



-M 



2 r sln0 <{) 



- 



r sin0 




f 


3^,M) ■ 


r sin0 


f 






e 

r sin0 



2 r sin0 0 ({)■ 
2 

1 + u, 

<v 






[Du" + ^ (3 - 2u")] 



= 0 



(H-16) 



the u terms are discarded and if the equation of 
is used, these equations become simply 



h'Q t 

3.(%) + — 0“r + 



2ip « 

3o(Du^) + ^DUg + 





(Du ) = 
r 


1 


1+T 


“ ¥ 


23„t 


1 + T 


II 

CD 

::J 

Q 


jL 

" ? 



e a^D 



(H-17) 
cont . 



- 9nD 
r 0 
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I 





(H-17) 



h p I ^ ^ ft 



1 .-M/2 1 . ^ 

TT 6 • — 7T 0 , D . 

4 r sin0 (j) 



In order to learn something about the behavior of the 
-M/2 h 

moment <e ^ sln9 vjhich appears in Eq. (H-13), 

consider the average of the third of Eq. (H-17) multiplied 
by h. 

P<h9oU(j)> = - ^ K'^> (H-18) 



r sln0 <}) 



(To obtain Eq. (H-l8), the higher order moments are 
neglec ..ed . ) . 

Because of the previous arguments concerning the rela- 
tionship of spatial and temporal partial derivatives, the 



following equation should hold for some f^ of order one. 



FTIm V/ " 



(H-19) 



If Eq. (H-19) is substituted into Eq. (H-13), the follow- 
ing equation results. 

-2 



<h9gd> = 



q * “Ifa 



p<hS n T> 



(H-20) 



If the changes of d and x occur over the same time, 
it must be true that 



<hd> = f p<hx> 

3 



(H-21) 



where f^ is a weak function of x® of order one (unless by 

chance f, + ^f is small). 

1 2 
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But in Eqs . (9-^) and (9-5), the correlations <td> 
and <OpT)d> are required. Since the ratio of D to D is 
(1 + u^/3 ) (see Eq. (9-6)) v.'hlch is very close to one, 

it can be assumed that for some functions and 3^.5 both 
vjeak functions of x° and of order one as described above, 
the follov/ing equations hold. 

<xd> = - 2 X^p<t^> 

<OoT)d>= - 3^pao<T^>. (H-22) 

If and 3^ are expanded in a Taylor series and only 
the constant terms are retained because of the supposed 
weak dependence of these functions on x°, the following 
approximate equations result. 

<Td> = - 2Xp<T^> 

<O..T)d>= - 3p3o<t"> . (H-23) 

These constants are then to be viewed as arbitrary parameters 
of the model. Eqs. (H-23) are therefore consistent with the 
conservation equations. 
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APPENDIX I 



The Dynamics of the Energy Momentum Tensor 
for a Matter Pilled Universe 

The dynamics of the energy-momentum distribution is 
given by 

T^^^; V = 0. (I-l) 



The conservation of energy equation 


corresponding to y 


equal to zero is given by Eq. (H-3) 


with u . and P equal to 
J 



zero. (It should be noted that Uj and P equal to zero with 
D unspecified are solutions to the three momentum equations 
for which y = 1, 2, 3- See Eqs . (H-l^l), (H-15), and (H-l6).) 



That is, 

3,D + |Oj,M)D= 0 


(1-2) 



3/2M 

If this equation is multiplied by e and integrated. 



the following equation results. 




^ g3M/2 


(1-3) 


With 

D = p + d 


(I-^) 


and 

^3/2M . ^3/20 (,+,)3/2 
equation (1-3) becomes for small t 


(1-5) 


p + d - ^3/2 g 2^^ ‘ 


(1-6) 
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Hence , 


P = ^ 3^/2 U-7) 


and 


d = - 1 pT . (1-8) 



That is, for a pressure free universe, the stochastic process 



d is a 


function of the stochastic prosess t. Hence, 

<Td> = - f P<T^> (1-9) 


and 


<Oj,x)d> = - 1 pd,<T^> . (I-IO) 
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APPENDIX J 



Solution of the Growth Equations 



for a Radiation Pilled Universe 



If y = <T^>, z = R** ^ , X = R/r'*’, the growth 



Eqs. (9-3^) and (9-35) are 




dx ^ 



(J-1) 



and 



^ + i z = ^ 

dx X 4 dx 

X 



(J-2) 



where R is a constant with the dimension of length. All 
the variables (y, z, x) in Eqs. (j-1) and (J-2) are 
dimensionless . 

The solution of these two coupled differential equations 
proceeds as follows. If Eq. (J-2) is multiplied by x**, 
the left hand side becomes a derivative of the function x“*z. 
This equation may then be Integrated directly to give 



From Eq. (J-3), the first and second derivatives of y 
with respect to x may be found in terms of z and x. If 
these are substituted into Eq. (J-1), a second order differ- 
ential equation in z results. ' 




1 



(J-3) 



g 




12^1 



2 d^Z 
dx ^ 



+ 9x 



dz 

dx 



32BXf^ 

+ 8(2 - X - 3)z = S. 



X 



(J-iO 



Equation (J--M can be transformed Into a second order 
differential equation v;lth constant coefficients by the 
change of variable 



X = e^ . (J-5) 

Thus , 

+ 8 ^ + 8(2 - X - 3)z = 323Xf e"^^ . (J-6) 

dt^ S 



The associated homogeneous equation of Eq. (j-6) is 

^-5. + 8 ^ + 8(2 - X - e)z = 0 . 



Its complementary solution z^ is 

z^ = Ae^^ + Be^^ (J-8) 

where i = - 4 h- 2 /2 ( X + 3 ) 
j = - 4 - 2/2(X+3) 



and X and 3 are arbitrary constants. 

The complementary solution z^ can be used to find the 
general solution to Eq. (J-6). This is done by assuming 
the general solution is given by Eq. (J-8) but v;ith A and 
B functiQHs of t. 

If this solution form is substituted into Eq. (J-6), 
a differential equation in A and b results. 



d-t 



2 ^ ^ 



+ 8(2 - X - 3) 



it 



(J-9) 

cont . 
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+ B [ 



(J-9) 



dt 



+ 8 ^ + 8(2 - A -B)] e 



J't 



r ^ d(e^^) ^ d(e^^) 

^ dt dt dt dt 



* at I? " i ] + 8[(e“) ff t 



= 32BXfg6 



-H 



(J-9) 



The first tv;o terms are zero by virtue of the fact that 
e^^ and e*^^ are solutions to Eq. (J-7). In the equation. 



z(t) = A(t)e^^ + B(t)e'^’^ 



(J-10) 



three dependent variables of t are present: z. A, and B. 

But in general, three variables can be made to satisfy three 
equations. Hov/ever, so far, only tvjo such equations have 
been specified, namely Eq. (J-^) and Eq . (J-10). Hence, 
a third equation can be freely chosen. A convenient choice 
is 



it ^ + gjt ^ ^ Q 



dt 



dt 



(J-11) 



so that Eq. (J-9) becomes 



d(e^^) ^ . d(e^^) OB 
dt dt dt dt 



32BAf„< 

g 



-H 



(J-12) 



These two equations can be written in matrix form, 



Zu = F 



(J-13) 



where Z = 



it 



Jt 



d(e^^) d(e^^) 



dt 



dt 
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u = 



dt 

dt 



P = 323Af e 

o 

By Cramer’s Rule, 



-i4t 



0 ' 



dt 



det Z 



det 
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dt 



det Z 



det 



32BXf„e 

g 



it 



d(e^^) 






so that 



dt 



6B 

dt 



32XBf^e' 

ti 



dt 



(i+^Ot 



J-i 



32X3^^© 



-(j+4)t 



J-i 



dCe-^’^) 

dt 



32X3fge 



-Ht 



(J-li^) 



(J-15) 



(J-16) 



for j7^i(X = -3) . 

Now, Eqs . (J-15) and (J-l6) may be integrated to give 



A = 






(j-i)(i+4) 



® (j-i)(J+T7 ® ^ 



(J--17) 



(J-18) 



where D and E are arbitrary constants. Hence, by Eq. (J-10) 



and by the identity x = e , 
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(J-19) 



2 = 



where F 

(J-3) , 

y 



- kXFx~^ + Dx^ + Ex-^* 

3f 

and is a new arbitrary constant. By Eq, 



„ ^ D i+4 ^ E j 

^ ^ T3 ^ ^ ^ 



+4 



(J-20) 



If these solutions for y and z are substituted into Eqs. 

(J-1) and (J-2), they are found to satisfy these equations. 

If y and z of Eqs. (J-19) and (J-20) are translated back 
, 2 

into <T > and <(9.t) >, Eqs. (J-19) and (J-20) become 



_ p , ^1/3 pi+^ , ^1/3 pJ+^ 

<T >.- Fi/ 3 + nrr ^ ^ ~W~ ^ 



^ - = - i]XF ,_R ^ + In + Dn /oR'^ 

a_ 1/3 1/3 1/3 



(J-21) 

(J-22) 



where i = - ^ + 2 / 2 ( X+3 ) 
j = - 4 - 2/2(X+3) 

and 1^/33 ^ 1 / 3 ’ ^ 1/3 arbitrary constants. ( 1^/3 ~ 

^/r+1+4 ^ E/R+'^‘^^.) 
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APPENDIX K 



Solution of the Grov\rth Equations for a Matter Filled Universe 

+ ^ 

If y = <T^>; z = — — <(3.t)^>; x = ^ (for some R"^), the 
growth Eqs. (9-65) and (9-66) are 



x^y " + |- xy ' 



. + 1 Z = 1 2_ 

* ^ x> 



3y = 2x^z 



(K-1) 

(K-2) 



If Eq. (K-1) is differentiated and Eq. (K-2) is substi- 
tuted into it, a third order differential equation for z 
results . 

2x^z’" + 27x^z" + 8^xz' + 42z = 0 . (K-3) 



If a change of variables is made such that 

X = e^ (K-^) 

then Eq . (K-3) becomes 

2z^^^ + 21z^^ + 6lz^ + il2z = 0 . (K-5) 

The general solution of this equation is simply, 

z = Ge“'^ + + Je"^^ . (K-6) 

In terms of x, Eq. (K-6) becomes 

z = Gx"^ + + Jx~^ . (K-7) 

If this equation is substituted into Eq. (K-2), an equation 
for y' results. If this equation is integrated, the 
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following expression for y appears. 




+ j Jx ^ + K 



(K-8) 



But Eqs. (K-7) and (K-8) are general solutions to the 
set of the tv/o coupled Eqs. (K-3) and (K-2). Pour integra- 
tions are required to solve this system v;hlch yields four 
arbitrary constants (G, H, J, and K above). However, the 
given system of equations (Eq. (K-1) and (K-2)) require only 
three integrations. Because Eq. (K-3) is derived from these 
equations by differentiation, Eqs. (K-7) and (K-8) must be 
the general solution to the given equations with one of the 
appropriate arbitrary constants equal to zero. By inspection, 
the equation 



although a solution to Eqs. (K-3) and (K-2), is not a solu- 
tion to the original equations. Hence, the general solution 
to Eqs. (K-1) and (K-2) must be 



y = K and z = 0 



1 — -7/2 — 

z = Gx + Hx '^ + Jx 



(K-9) 




(K-10) 



Hence 




(K-11) 



and 




(K-12) 
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APPENDIX L 



Classification of Solutions for the Radiation Filled Universe 



2 

With the exception of the integration of R' , Eqs . (E-20) 
and (E-21) of Appendix E give the solutions to Eqs. (E-1) 
and (E-2). These solutions are, with k neglected as usual 
(see Eqs. (9-^1) to (9-^5)), 



where 



Kp = - 2ir. f r">A 
R* 



R'^ =aj,R"^ + 



0(1, 73) r 



h/3 



lOX + - 1 + 



C(Di/3) 



■ D^/3 12 



lOX + - 1 - 



■1/3 



Cdi/3) ^ 

—3^ 



2X + /2(X+B) 
/2(X+B) 



"1/3 " 



C(Di/ 3 ) j- 2 x _ d(X+3) 



/2(X+B) 



m. 



■1/3 



= i = - + 2/2 (X + B) 



= j = - - 2/2 (X + B) 

^1/3, 



2B ■ 

/2(X+B) 

2B ■ 

/2(X+B) . 



(L-1) 

(L-2) 

(L-3) 

iL-h) 

(L-5) 

(L-6) 

(L-7) 

(L-8) 



The A subscript designates the type of fluctuation. (I2/3 ~ 
type corresponds to a situation v;here the nonuniformity 
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type corresponds to 



Increases v^lth increasing R and ~ 

a situation where the nonuniformity decreases v/lth Inceasing 

R. ) 

There are four kinds of solutions based upon the manner 
in which <p and R are displaced from the fluctuation free 
solutions Kp* and R* . That is, f^ and can in principle 
be each positive or negative depending on the choice of X 
and 3 . 

In order to determine how a given (A, 3) displace Kp and 
R from Kp* and R* , it is necessary to locate those values 
(X,3) for v.’hich f and b are equal to zero. Then a set (A, 3) 
which lies on a particular side of the two curves of zero 
f^ and zero b^ determines the sign of f^ and b^ and hence 
the direction of displacement of <p and R from Kp* and R* . 

In order to find the line of zero f^ in A-3 space, con- 
sider first the following definition. 

n = /2(A+3) . (L-9) 

The set (A, 3) for v;hich f^ equals zero is given by the 
equation 

lOA + ^3 - 1 = + - -- (L-10) 

/2(A+3) 

where the plus sign is used if A represents 

minus sign if A represents This equation as it stands 

is rather difficult to plot directly. But the plotting is 
simplified if Eq. (L-10) is separated into two parametric 
equations v.^ith parameter n. 
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J 









The two parametric equations 3 (n) and A(n) are found to 



be 



p _ n [ 5n^ - 1 
^ 2 [ 3n + 1 



(L-11) 



n 2n^ + n - 1 

2 3n + 1 



(L- 12 ) 



where n > 0 . The top sign corresponds to the ~ type 

fluctuation and the bottom sign to the I2/3 ~ type 
fluctuation . 

By varying n, the equation 3 (A) can be traced out. This 
equation is plotted in Figs. 3 and ^ for both Ij/3 ^ 1/3 
type fluctuations respectively. 

The equation of A vs . 3 Tor b^ = 0 is much simpler to 
find. With straightforward algebra, this equation is 



3 = A( 2 A - 1 ) (L-13) 

for both the 1^/3 ^ 1/3 ~ fluctuations. This equa- 

tion is also plotted in Figs. 3 and . It should be noted 
that Eq. (L-I3) is quadratic in A. Hence, according to this 
equation, there should be two values of A corresponding to 
a given value of g. Hov;ever, only one value satisfies the 
requirement that b^ equal zero. That there are tv/o values of 
A comes from the fact that Eq. (L-I3) is derived by a squar- 
ing operation v/hich thereby introduces an extra root. 
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APPENDIX M 



Index and Location of Symbols Used in Dissertation 

The follov;ing list contains the symbols used in this 
work. Accompanying each symbol is the location v/here the 
definition or explanation of the symbol may be found. 



1 . 


L*, G, c 




I 


2 . 


ds, g^^, x^, X®, r, 0 , d. 


da^, Rq, 


IV-A 




Rj^, k, e°, e^, M 




IV-A 


3. 






IV- B 


ll. 


tpv,3J,3^, C- V- 

pX 

G j K 

V ’ 


R^ , R, 

V ’ ’ 


V-A 


5. 


D, P, u, D, P, a 




V-F 


6 . 


R 




VI-A 


7. 


Pj d 




IX-A, B, 


8 . 


M, V, AD, am, AV, g, V, X 


, 8 


IX-A-2 


9. 


^1/3’ ^1/3’ ^1/3’ 




IX-A-3 


10 . 


Cd^/ 3 ), CCD^/^), C(A) 




IX-A-k-a 


11 . 


p*, R* 




IX-A-4-b 


12 . 






IX-B-3 


13. 


B(IJ, B(DJ 




IX-B-4 


IH. 


6p, D, I, Tj, Tg, <xd.. 




X-A 



13 ^ 



15. 



X-B 



<5T, T, Hp, L, 6p , p , c, R„, R , R 
’ R’ ’ ^o’^o’ ’ o’ max’ e 



16 , 

17 

18 

19 

20 
21 



V, V, V, P, p 

P(6;x^), f(D;x^), h[D(x^)] 

^ V ’ ’ ^v’ ' 



P^, p^, cx^, V, in, K, ('), C, B, a, a^ , 
a^, f, b 



P> b., h, Pj> ^ 2 > ^3’ 



y, X, z, R , t 



22. t, A, B, Z, u, F, Pg> ^5 P 

23. G, H, J, K 



2^ • > Pf- > Pj-\ > by , bj.. , m_ , 

±1/3 -^ 1/3 ^ 1/3 -^ 1/3 



A 

B 

D 

E 

H 

J, K 

J 

K 



% ’ ^A’ ^A’ "" 

Ui/3 A A 
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Figure 1. Classification Designations 
for the Increasing Type Fluctuations. 

The solid lines are the fluctuation free solutions for R 
and p. The broken lines depict the possible deviations of 
R and p from these unperturbed solutions. The letters K, 

L, M, and N label the four possible combinations of perturbed 
solutions for R and p from their unperturbed values . 
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Figure 2. Classification Designations 
for the Decreasing Type Fluctuations. 

The solid lines are the fluctuation free solutions for R 
and p. The broken lines depict the possible deviations of 
R and p from these unperturbed solutions. The letters K, L, 
Mj and N label the four possible combinations of perturbed 
solutions for R and p from their unperturbed values. 
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Figure 3* Solution Types in X-3 Space for the Increasing 
Type Fluctuations for the Radiation Filled Universe 
The plane is sectioned into four regions. The perturbed 

splutions for each region for R and p are depicted as broken 
lines. The fluctuation free solutions are represented as 
solid lines. 
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Figure Solution Types In X-$ Space for the Decreasing 

Type Fluctuations for the Radiation Filled Universe. 

The X-3 plane Is sectioned Into three regions. The perturbed 
solutions for each region for R and p are depicted as broken 
lines. The fluctuation free solutions are represented as 
solid lines. 
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Figure 5. Lines of Constant i and j. 
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